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ABSTRACT 

We present a simple, largely empirical but physically motivated model to interpret 
the mid- and far-infrared spectral energy distributions of galaxies consistently with 
the emission at ultraviolet, optical and near-infrared wavelengths. Our model relies on 
an existing angle-averaged prescription to compute the absorption of starlight by dust 
in stellar birth clouds and in the ambient interstellar medium (ISM) in galaxies. We 
compute the spectral energy distribution of the power reradiated by dust in stellar 
birth clouds as the sum of three components: a component of polycyclic aromatic 
hydrocarbons (PAHs); a mid-infrared continuum characterising the emission from hot 
grains at temperatures in the range 130-250 K; and a component of grains in thermal 
equilibrium with adjustable temperature in the range 30-60 K. In the ambient ISM, 
we fix for simplicity the relative proportions of these three components to reproduce 
the spectral shape of diffuse cirrus emission in the Milky Way, and we include a 
component of cold grains in thermal equilibrium with adjustable temperature in the 
range 15-25 K. Our model is both simple and versatile enough that it can be used to 
derive statistical constraints on the star formation histories and dust contents of large 
samples of galaxies using a wide range of ultraviolet, optical and infrared observations. 
We illustrate this by deriving median-likelihood estimates of the star formation rates, 
stellar masses, effective dust optical depths, dust masses, and relative strengths of 
different dust components of 66 well-studied nearby star-forming galaxies from the 
Spitzer Infrared Nearby Galaxy Survey (SINGS). We explore how the constraints 
derived in this way depend on the available spectral information. From our analysis 
of the SINGS sample, we conclude that the mid- and far-infrared colours of galaxies 
correlate strongly with the specific star formation rate, as well as with other galaxy- 
wide quantities connected to this parameter, such as the ratio of infrared luminosity 
between stellar birth clouds and the ambient ISM, the contributions by PAHs and 
grains in thermal equilibrium to the total infrared emission, and the ratio of dust mass 
to stellar mass. Our model can be straightforwardly applied to interpret ultraviolet, 
optical and infrared spectral energy distributions from any galaxy sample. 
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1 INTRODUCTION 

The spectral energy distributions of galaxies contain valu- 
able information about their contents in stars, gas and 
dust. Direct ultraviolet, optical and near-infrared radiation 
from stars provides clues on the past star formation his- 
tory, chemical enrichment and attenuation by dust. Neb- 
ular emission lines produced by the gas heated by young 
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stars provide further clues on the current star formation ac- 
tivity and the physical conditions of the star-forming gas. 
At wavelengths A > 3 /tm, the mid- and far-infrared emis- 
sion reflects the heating of dust in different components 
of the interstellar medium (ISM) by stars of all ages. Ob- 
servations at ultraviolet, optical and infrared wavelengths 
are now becoming available for large samples of galax- 
ies. These include data collected in the ultraviolet by the 
Galaxy Evolution Explorer ( GALEX, \Martin et al.ll2005h . in 
the optical by the Two-degree Field Galaxy Redshift Sur- 
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vey llColless et al. 2001) and the Sloan Digital Sky Sur- 
vey (SPSS. Istoughton et al.ll2002l ). in the near-infrared by 
the T wo Micron All Sky Survey (2MASS, ISkrutskie et all 
1 19971 ). in the mid- and fa r-infrared by the Inf rared As- 
tronomical Satellite (IRAS, iBeichman et a.1. Il988h . the In- 
frared Space Observator y (ISO, Kessler et al.l l l996l ) and the 
Spitzer Space Telescope l| Werner et al.ll2004 ). and in the sub- 
millimetre by the Sub-millimeter Common User Bolome- 
ter Ar ray (SCUBA) on t he James Clerk Maxwell Tele- 
scope l|Holland et al.lll999l ). Extracting constraints on the 
stellar populations and ISM of galaxies from these multi- 
wavelength observations requires the consistent modelling 
of the emission by stars, gas and dust. 

A standard approach to model consistently the emission 
from stars and dust in galaxies has been to solve the radia- 
tive transfer equation for idealised (bulge + disc) spatial 
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distributions of stars and dust 
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Misselt et all l200ll ; iPopescu et al.1 l2000l : iMisiriotis et all 
200 ll ). Early models of thi s type did not inc lude the evolu- 



tion of stellar populations. ISilva et al] (|l998l ) were the first 
to couple radiative transfer through a dusty ISM and the 
spectral (and even chemical) evolution of stellar populations. 
Their model also accounts for the fact that stars are born 
in dense molecular clouds, which dissipate after some time, 
and hence, that newly bor n stars are more attenuated than 
older stars (see also, e.g., ICharlot fc Fallll2000l ; iTuffs et all 
120041 ). This type of sophisticated model is useful to interpret 
in detail the emission from individual galaxies in terms of 
constraints on stellar populations and the spatial distribu- 
tion and physical properties of the dust. However, because of 
the complexity of radiative transfer computations, it is not 
optimised to derive statistical constraints from observations 
of large samples of galaxies. 

A more rece nt mo del of starburst galaxies by 
iDopita et al.l l|2005l . see also lGroves et alj|2007l ) incorporates 
the consistent treatment of the spectral evolution of stellar 
populations, the dynamic expansion of H 11 regions and ra- 
diative transfer of starlight through gas and dust. The au- 
thors of this model provide a simple parameterization of the 
ultraviolet, optical and infrared spectra of starburst galaxies 
by adding the spectra of different types of compact H 11 re- 
gions and their surrounding photo-dissociation regions. This 
model provides a fast and flexible tool to interpret starburst 
galaxy spectra in terms of the physical parameters of star- 
forming regions. However, it is not designed to be applicable 
to more quiescent galaxies, in which older stellar populations 
dominate the emission. 

In parallel to these theoretical studies, observations by 
IRAS and ISO have motivated the development of sim- 
ple, empirically calibrated spectral libraries to interpret 
the infrared emission from gala xies at waveleng t hs bet ween 
3 and 1000 fim. For example, ICharv fc Eibazl (|200ll ) and 
iDale fc Heloul (|2002l ) both present single-parameter fami- 
lies of infrared spectra to relate an observed spectral en- 
ergy distribution to either the total infrared luminosity 
of a galaxy or the intensity of the interstellar radiation 
field heating the dust. These libraries can be applied eas- 
ily to the interpretation of large galaxy samples. They 
have proved useful to characterise the infrared emission 
from observed galaxies and to inv estigate the origin of the 
cosmic infrared background (e.g. iFranceschini et all l200ll ; 



20041 ; IDale et al.ll2005l : iMarcillac et aij l2006h . A disadvan- 
tage of this approach is that it does not relate consistently 
the infrared emission of the dust to the emission from stel- 
lar populations. Another potential limitation is that most 
existing spectral libraries were calibrated using local galaxy 
samples, and hence, they may not be applicable to studies 
of the infrared emission from galax ies at all redshifts (e.g. 
iPope et al.ll2006l ; IZheng et al.1l2007l ). 

In this paper, we present a simple, largely empirical 
but physically motivated model to interpret the mid- and 
far-infrared spectral energy distributions of galaxies consis- 
tently with the emission at ultraviolet, optical and near- 
infrared wavelengths. We comp ute the spectral ev o lution 
of stellar populations using the iBruzual fc Charlotl (|2003h 
population synthesis code. To describe the attenuation of 
starlight by du s t, we appeal to the two-component model of 
ICharlot fc Fail l|2000h . This has been shown to account for 
the observed relations between the ultraviolet and optical 
(line and continuum) spectra and the total infrared lumi- 
nosities of galaxies in wide ranges of star-formation activ- 
ity a nd dust content (|Brinchmann et all 120041 : iKong et al.l 
2004). We use this model to compute the luminosity ab- 
sorbed and re-emitted by dust in stellar birth clouds (i.e. gi- 
ant molecular clouds) and in the ambient (i.e. diffuse) ISM in 
galaxies. We then distribute this luminosity in wavelength to 
compute infrared spectral energy distributions. We describe 
the infrared emission from stellar birth clouds as the sum 
of three components: a component of polycyclic aromatic 
hydrocarbons (PAHs); a mid-infrared continuum character- 
ising the emission from hot grains at temperatures in the 
range 130-250 K; and a component of grains in thermal equi- 
librium with adjustable temperature in the range 30-60 K. 
In the ambient ISM, we fix for simplicity the relative pro- 
portions of these three components to reproduce the spectral 
shape of diffuse cirrus emission in the Milky Way, and we 
include a component of cold grains in thermal equilibrium 
with adjustable temperature in the range 15-25 K. 

This simple but versatile model allows us to derive sta- 
tistical estimates of physical parameters such as star for- 
mation rate, stellar mass, dust content and dust properties, 
from combined ultraviolet, optical and infrared galaxy spec- 
tra. To achieve this, we adopt a Bayesia n approach simi- 
lar to t hat successfully employed by, e.g., Kauffmann et al.l 
ll2003h.lBrinchmann etafl (|2004l ), lGallazzi et all i|2005l ) and 
ISalim et al.l ( 20071 ) to interpret ultraviol et, optical and near- 
infrar ed ga laxy spectra us i ng on ly the IBruzual fc Charlotl 
l|2003h and ICharlot fc Fail (|2000l ) models. As an example, 
we derive median-likelihood estimates of a set of physi- 
cal parameters describing the stellar and dust contents of 
66 star-forming galaxies from the Spitzer Infr ared Nearby 
Galaxy Survey (SINGS. iKennicutt et aill2003l ). Our model 
reproduces well the observed spectral energy distributions 
of these galaxies across the entire wavelength range from 
the far-ultraviolet to the far-infrared, and the star forma- 
tion histories and dust contents of the galaxies are well con- 
strained. We explore how the constraints derived in this way 
depend on the available spectral information. From our anal- 
ysis of the SINGS sample, we conclude that the mid- and 
far-infrared colours of galaxies are tightly related to the spe- 
cific star formation rate and to other galaxy-wide properties 
connected to this parameter. 
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We present our model of the combined ultraviolet, op- 
tical and infrared spectral energy distributions of galaxies 
in Section [2] In Section [3] we first describe our approach to 
derive statistical constraints on galaxy physical parameters 
from multi-wavelength observations. Then, we use this ap- 
proach to interpret observations of the SINGS galaxy sam- 
ple taken with GALEX, 2MASS, Spitzer, ISO, IRAS and 
SCUBA. We compare our results to those that would be ob- 
tained using previous prescriptions of the infrared emission 
from galaxies. We also discuss potential sources of system- 
atic errors. Section [4] summarises our conclusions. 



2 THE MODEL 

In this section, we describe our model to compute the mid- 
and far-infrared spectral energy distributions of galaxies 
consistently with the emission at ultraviolet, optical and 
near-infrared wavelengths. In Section 12.11 we first briefly 
review the stellar population synthesis code and the two- 
component dust model used to compute the emission of 
starlight and its transfer through the ISM in galaxies. Then, 
in Section 12.21 we present our model to compute the spec- 
tral energy distribution of the infrared emission from dust. 
We calibrate this model using a sample of 107 nearby star- 
forming galaxies observed by IRAS and ISO. In Section [2. 31 
we show examples of combined ultraviolet, optical and in- 
frared spectral energy distributions of different types of star- 
forming galaxies. 



2.1 Stellar emission and attenuation by dust 

We use the latest version of the lBruzual fc Chariot! (|2003l ) 
stellar population synthesis code to compute the light pro- 
duced by stars in galaxies. This code predicts the spectral 
evolution of stellar populations at wavelengths from 91 A to 
160 /im and at ages between 1 x 10 5 and 2 x 10 10 years, 
for different metallicities, initial mass functions (IMFs) 
and star formation histories. We use the most recent ver- 
sion of the code, which incorporates a new prescription by 
iMarigo fc Girardil (120071 ) for the TP-AGB evolution of low- 
and intermediate-mass stars. The main effect of this pre- 
scription is to improve the predicted ne ar-infrared col ours 
of intermediate-age stellar populations (jBruzuall 120071 ; see 
also Chariot & Bruzual, in preparation). In all applications 
throughout this paper, we adopt the Galactic-disc IMF of 
IChabrierl l|2003l ). 

We compute the attenuation of starlight by du s t usin g 
the simple, angle-averag ed model of lCharlot fc Fall! (2000). 
This accounts for the fact that stars are born in dense 
molecular clouds, which dissipate typically on a timescale 
to ~ 10 7 yr. Thus the emission from stars younger than 
to is more attenuated than that from older stars. Following 
ICharlot fc Fall (|200Ch . we express the luminosity per unit 
wavelength emerging at time t from a galaxy as 



(1) 



Jo 

where ip(t — t') is the star formation rate at time t — t', S\(t') 
is the luminosity per unit wavelength per unit mass emitted 
by a stellar generation of age t' , and T\(t') is the 'effective' 
absorption optical depth of the dust seen by stars of age t' 



(i.e. averaged over photons emitted in all directions by stars 
in all locations within the galaxy). The time dependence of 
T\ reflects the different attenuation affecting young and old 
stars in galaxies, 



for t' < t , 
for t' > t . 



(2) 



{~ BC , -ISM 

Here t a bc is the effective absorption optical depth of the dust 
in stellar birth clouds and f A t hat in the amb i ent IS M. 
We also adopt the prescription of ICharlot fc Falll (|2000i ) to 
compute the emergent luminosities Lnaif) and L-g_p(€) of the 
Ha (A = 6563 A) and H/3 (A = 4861 A) Balmer lines of hy- 
drogen produced by stars in the birth clouds. This assumes 
case B recombination and includes the possible absorption 
of ionising photons by dust before they ionise hydrogen. 

The shape of the effective absorption curve depends on 
the combination of the optical properties and spatial distri- 
bution of the dust. We adopt the following dependence of 
f BC and f A ISM on wavelength: 



t x BC = (1 - v)tv (A/5500. 
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(3) 
(4) 



where fv is the total effective V-band absorption optical 
depth of the dust seen by young stars inside birth clouds, and 
is the fraction of this contributed by 



/'• : 



dust in the ambient ISM. The dependence of f A ISM on A" 7 
is well constrained by the observed relation between ratio of 
far-infrared to ultraviolet luminosity and u ltraviolet spec- 
tral s lope for nearby starburst galaxies fsee lCharlot fc Falll 
2000). The dependence of r x BC on wavelength is less con- 
strained by these observations, because stellar birth clouds 
tend to be optically thick, and hence, stars in these clouds 
contribute very little to the emerge nt radiation (ex c ept in 
the emission lines). For simplicity, ICharlot fc Fan] (2000) 
adopt f A BC ex A -0 ' 7 by analogy with f A ISM . We adopt here 
a slightly steeper dependence, f x BC oc A -1,3 (equation [3]), 
which corresponds to the middle range of the optical prop- 
erties of dust grains between the Milky Way , the Large and 
the S mall Magellanic Clouds (see section 4 of lCharlot fc Fall 
2000). This choice is motivated by the fact that giant molec- 
ular clouds can be assimilated to foreground shells when at- 
tenuating the light from newly born stars. In this case, the 
effective absorption curve should reflect the actual optical 
properties of dust grains. We emphasise that the depen- 
dence of r BC on wavelength has a negligible influence on 
the emergent ultraviolet and optical continuum radiation. 
It affects mainly the attenuation of emission lines in galax- 
ies with large t bc /t^ m and h ence small /x (Section 13.4.21 
below; see also lWild et~ai1l2007f l. 

The fraction of stellar radiation absorbed by dust in the 
stellar birth clouds and in the ambient ISM is reradiated in 
the infrared. We write the total luminosity absorbed and 
reradiated by dust as the sum 

Lrw=id c w+irw, (5) 



where 



(6) 



d\(l-e~ f * ) J dt'ip(t-t')Sx(t') 

is the total infrared luminosity contributed by dust in the 
birth clouds, and 
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T ISM/ 
^d 



poo . pt 

(t) = / dA(l-c^ SM / dt'ip(t-t')Sx(t') 



(7) 



is the total infrared luminosity contributed by dust in the 
ambient ISM. For some purposes, it is also convenient to de- 
fine the fraction of the total infrared luminosity contributed 
by the dust in the ambient ISM: 



u(t)=Li sM ( t )/Lr(t). 



(8) 



This depends on the total effective V-band absorption opti- 
cal depth of the dust, fv, the fraction p of this contributed 
by dust in the ambient ISM, and the star formation history 
(and IMF) determining the relative proportion of young and 
old stars in the galaxy. 



2.2 Infrared emission of the dust 

We now present a simple but physically motivated prescrip- 
tion to compute the spectral distribution of the energy rera- 
diated by dust in the infrared (i.e., the distribution in wave- 
length of Lj 30 and Lj SM ). By construction, the infrared 
emission computed in this way can be related to the emission 
at shorter wavelengths using equations (HJ — (EJ above. 



2. 2. 1 Components of infrared emission 

The infrared emission from galaxies is generally attributed 
to three main constituents of interstellar dust: polycyclic 
aromatic hydrocarbons (PAHs) , which produce strong emis- 
sion features at wavelengths between 3 and 20 /im; 'very 
small grains' (with sizes typically less than 0.01 /im), which 
are stochastically heated to high temperatures by the ab- 
sorption of single ultraviolet photons; and 'big grains' (with 
sizes typically between 0.01 and 0.25 /im), which are in ther- 
mal equilibrium with the radiation field. This picture arises 
from detailed models of the sizes and optical properties of 

3 ISM of the Milky 
_ Mathis et al.1 1 19771 : 
iDraine fc Leell 19841 : iLeeer fc Pueetlll984l ). Here, we build on 
the results from these studies to describe the different com- 
ponents of infrared emission in galaxies, without modelling 
in detail the physical properties of dust grains. As men- 
tioned in Section [1] the motivation for this approach is to 
build a model simple enough to derive statistical constraints 
on the star formation and dust properties of large samples 
of observed galaxies, based on consistent fits of the ultravio- 
let, optical and infrared emission. We describe the different 
components of infrared emission in our model as follows: 



dust grains and PAH molecules in th e ISM of the Milky 
Way and other nearby galaxies (e.g., 



(i) PAH and near-infrared continuum emission. The 
mid-infrared spectra of most normal star-forming galax- 
ies are dominated by strong emission features at 3.3, 6.2, 
7.7, 8.6, 11.3 and 12.7 /im. Although still uncertain, the 
carriers of these features are generally accepted to be 
PAH molecules transiently excited to high internal energy 
levels after the absorption of single ultraviolet photons 
ilLeger fc PugetJll984l: lAllamandola et al.lll985l : iLeger et ail 
ll989l : lAllamandola et all 19991 ). PAH emission tends to peak 
in the 'photo-dissociation regions' (PDRs) at the inter- 
face between ionised and molecular gas in the outskirts of 
Hn regions, where PAH molecules c an survive and tran- 
sient heating is most efficient (e.g., ICesarskv et al.l 1 19961 : 



IVerstraete et al.lll996l : iRapacioli et al.| [2005). In these envi- 
ronments, the non-ionising ultraviolet radiation from young 
stars dominates the energy balanc e and can dissociate 
mole cules such as H2 and CO (see iHollenbach fc Tielend 
Il997l for a review) . 

Observations with /SO/ISOPHOT of the mid-infrared 
spectra of 45 nearby, normal star-forming galaxies do not 
reveal any strong variation of the PAH emission spectrum 
with galaxy properties su ch as infrared colo urs and infrared- 
to-blue luminosity ratio (He lou et alJ liToOO), For simplicity, 
we adopt a fixed template spectrum to describe PAH emis- 
sion in our modelQ We use the mid-infrared spectrum of the 
photo-dissociation region in the prototypical Galactic star- 
forming region Ml 7 SW extract e d from the JffO/ISOCAM 
observations of ICesarskv et~aH (I1996T ) by iMadden et al.l 
(2006). The observed spectrum does not extend blueward 
of 5 /im. We extend it to include the 3.3 /im PAH emis- 
sion feature using the Lorentzian profile parametrised by 
IVerstraete etafl feoOll ). F^ orentz = / [1 + (x - xof/a 2 ]' 1 , 
where x — 1/A is the wavenumber, xo = 3039.1 cm -1 the 
central wavenumber of the feature and a — 19.4 cm - half 
the FWHM. We set the amplitude fo so that the luminos- 
ity of the 3.3 /im emission feature is 10 per ce nt of that of 
the 1 1.3 fim feature (the relative ratio found bv lLi fc Drainel 
l200ll for neutral PAHs). We write the spectral energy dis- 
tribution of PAHs in our model 



,PAH _ r M17 



d\ L\ 



(9) 



where 17 is the adopted M17 spectral template, and /^ AH 
is normalised to unit total energy. 

The above /SO/ISOPHOT observations of a sample of 
45 normal star-forming galaxies also reveal a component of 
near-infrared continuum emission characterised by a very 
high colour tempera ture (~ 1000 K) at wavelengths be- 
tween 3 and 5 fim dLu et alJliool ). This component ac- 
counts typically for at most a few percent of the total in- 
frared luminosity, but it contributes significantly to the ob- 
served Spitzer/IRAC band fluxes at 3.6 and 4. 5 /im. It is also 
prese nt in the spectra of reflection nebulae (|Sellgren et al.l 
Il983h and in the d iffuse cirrus emissio n of the Milky Way 
l|Dwek et ail Il997l ; iFlagev et alll2006h . The origin of this 
emission is still uncertain, and it could be related to the 
stochastic heating o f PAH m o lecule s or carbon grains (e.g. 
IFlagev et al]|200d ). iLu et all (|2003T ) find that the strength 
of this emission correlates well with that of PAH features in 
the spectra of the galaxies in their sample. In particular, the 
ratio of the continuum flux density at 4 fim to the mean flux 
density of the 7.7 fim PAH feature has a roughly constant 
value of 0.11. 

To implement this component of near-infrared continuum 
emission associated to PAH emission in our model, we use 
a gre ybody (i.e. modifi ed blackbody) function of the form 



a grey body (i.e. modin 
fe.g.. lmidebrandlll983h 



1 In reality, changes in PAH molecule size and ionisation state, 
variations in metallicity and contamination by an AGN could 
cause variations of up to 40 pe r cent in the rela t ive strengths 
of some PAH emission features llHelou et alj|2000l : iBrandl et alj 
120061: ISmith et al1l2007l; iGalliano et al.ll2007ft . We ignore this re- 
finement here. 
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k x B\(T d ) 



d\n x B x (T d ) 



(10) 



where B\(Td) is the Planck function of temperature Td, K\ 
is the dust mass absorption coefficient, and l^ d is normalised 
to unit total energy. The dust mass absorption coefficient is 
usually approximated by a single power law of the form 



k\ ex A 13 . 



(11) 



where /3 is the dust emissivity index. Models and observa- 
tions of infrared to sub-millimetre dust emission and labora- 
tory studies of carbonaceous and amorphous silicate grains 
suggest values of f3 in the range 1 < f3 < 2 , with some de- 
pendence on grai n size and temperatu r e (lAndriessel 



Hildebrandl 1 19831; IPraine fc Led 1 1984 iReach et al 



Aglad ze et alj 1 19961 ; iBoulanger et al 



1996 



19981; see also section 3 4 of lDunne fc Ealed,2001, and section 



1974 



1995 



Mennella et alj 



Dale fc Heloul [2002 for discussions on the constraints i 
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/?). Typically, these studies favour /3 ~ 1 for small carbona- 
ceous grains, which radiate most of their energy at mid- 
infrared wavelengths, and /3 « 1.5 — 2 for big silicate grains, 
which reach lower temperatures and radiate most of their 
energy at far-infrared and sub-millimeter wavelengths. We 
therefore adopt = 1 in equation (|lip to compute the near- 
infrared continuum emission in our model. We scale this 
emission so that the continuum flux density at 4 fim be 0.11 
times the mean flux density of the 7.7 /jm PAH feature (see 
above). We find that a temperature Tj = 850 K provides 
optimal fits to the observed spectral energy distributions 
of galaxies at wavelengths between 3 and 5 /im (see, e.g., 

Fig.ru. 

(ii) Mid-infrared continuum emission from hot dust. In 
addition to PAH features, the mid-infrared spectra of star- 
forming galaxies (out to wavelengths A ~ 40 ^m) also in- 
clude a component of smooth continuum emission. This 
component is generally attributed to a continuous distri- 
bution of small grains with very low heat capacity, which 
are stochastically heated to high temperat ures by the ab - 
sorption of single ultraviolet photons (e.g., Sel lgrenlll984h . 
The accurate modelling of this emission would require the 
computation of temperature fluctuations of grains of differ- 
ent sizes and optical properties caused by the absorption of 
ultr aviolet photons in different interstellar r adiation fields 
(e.g. |Purcell|[l976l ; lAannestad fc Kenvonlll979l ) . For simplic- 
ity, in our model, we describe the 'hot' mid-infrared contin- 
uum emission as the sum of two greybodies (equation [TU| of 
temperatures Td = 130 and 250 K, with equal contributions 
to the total infrared luminosity, 



7 MIR _ /.250 | ;130\ 
'A — (l\ + 'A ) 



dA(Zf + O 



(12) 



We find that these two temperatures reproduce in an opti- 
mal way the range of infrared colours of star-forming galax- 
ies (see Section T2 . 2 . 2 1 below) . 

(iii) Emission from grains in thermal equilibrium. At far- 
infrared wavelengths, the emission from galaxies is gener- 
ally dominated by dust grains in thermal equilibrium at low 
temperatures. The grain temperatures depend sensitively on 
the intensity of the interstellar radiation field. This is why 
the peak of the far-infrared spectral energy distribution of 
a galaxy is a good indicator of dust heating in the ISM. 

We consider two types of grains in thermal equilibrium 



in our model: warm grains, which can reside both in stel- 
lar birth clouds and in the ambient ISM, with characteristic 
temperatures T^ c and T^ M , respectively; and cold grains, 
which can reside only in the ambient ISM, with characteris- 
tic temperature Tq Sm . To account for the observed disper- 
sion in the peak of the far-infrared emission of galaxies with 
different star formation activities, we allow the warm-grain 
temperature to vary between 30 and 60 K and the cold-grain 
temperature to vary between 15 and 25 K. We describe the 
emission from grains in thermal equilibrium using greybody 
spectra fequation llO[) with emissivity index (5 — 1.5 for warm 
dust and (3 — 2.0 for cold dust (equation I X 1 p . 

In summary, the infrared spectral energy distribution 
of stellar birth clouds in our model can be written 

7-BC_/eBC ; PAH f BC , MIR , ( BC ; T W°W 1 f \ T tot f-|o\ 
L A,d - (?PAH*A +?MIR'A +?W 'a ) — Jt*) L 'd A 16 ) 

where L d ot is the total infrared luminosity reradiated by 
dust (equation [5} , / M is the fraction of this contributed by 

rp BC 

the ambient ISM (equation [Sj , ^ AH , 1™ IR and l x w are 
computed using equations (|9}, (|12|l and (|10[) . and £pah, £mir 
and are the relative contributions by PAHs, the hot 
mid-infrared continuum and grains in thermal equilibrium 
to the total infrared luminosity of the birth clouds. These 
satisfy the condition 



t BC t BC t BC 
SPAH T ?MIR T ?W 



1 . 



(14) 



Similarly, the infrared spectral energy distribution of 
the ambient ISM can be written 



r ISM 

-^A,d 



_ fVISM / PAH , C ISM / MIR , C ISM ] L 
— (SPAH 'A T SMIR 'A T ?W 'a 



rr ISM 
+ th SM ll C ) 



r r tot 



(15) 



where L 



I 



and l x c 



are computed using 



(16) 



equations ©, CE1 and dlOj) , and £pah, Cmir, £w and £ C SM 
are the relative contributions by PAHs, the hot mid-infrared 
continuum and warm and cold grains in thermal equilibrium 
to the total infrared luminosity of the ISM. These satisfy the 
condition 

C ISM _i_ t ISM , C ISM , C ISM _ i 
?PAH ' ?MIR T ?W T ?C — 1 • 

In practice, we can fix the shape of the mid-infrared 
spectral energy distribution of the ambient ISM in our model 
to keep the number of adjustable parameters as small as pos- 
sible. This is justified by the fact that the intensity of the 
average radiation field heating dust in t he diffuse ISM of nor- 
mal galaxies is roughly constant (e.g., Hclou 1 9861). More- 
over, sophi sticated mod e ls of dust emission bv lLi fc Draind 
(|200lh and lDraine fe~Ll (|2007i ) suggest that even large vari- 
ations of the intensity of the interstellar radiation field have 
only a small influence on the overall shape of the diffuse mid- 
infrared spectral energy distribution. In these conditions, 
we can appeal for example to observations of high-Galactic- 
latitude (cirrus) dust emission in the Milky Way to constrain 
the mid-infrared spectral energy distribution of the ambient 
ISM in our model. In Fig. [T] we show the model spectral 
energy distribution L\ ^ (in black) computed using equa- 
tion (|15p that best fits the mean Galactic cirrus emission ob- 
served by the Cosmic Background Explorer/Diffuse Infrared 
Background Experiment ( C0-B-5/DIRBE) at wavelengths 
between 3.5 and 240 /im l|Dwek et al.lll997l ). The blue lines 
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Figure 1. Best model fit (in black) to the observed mean spectral energy distribution of the Galactic cirrus emission. The model was 
co mputed using equa tion 11511 and the parameters listed in equation dlTI . The red fill ed squares ar e the COBE/DIRBE observations 
of iDwek et al.l dl997T ). Also sho wn for refer e nce a re the Spitzer /IRAC observations of iFlagev et al.l ||2006| . red open squares) and the 
COBE/FIRAS observations of IDwek etall Jl997l . green line). The blue lines show the decomposition of the model in its different 
components (see Section 12.21 1: near-infrared continuum (dotted); PAHs (short-dashed); hot mid-infrared continuum (long- dashed) and 



warm and cold grains in thermal equilibrium (solid). The fit residuals (V 
model fluxes in a given photometric band, are shown at the bottom 



obs 



r mod 
A 



where £^ bs and L™ od are the observed and 



show the different components of the best-fit model, with 
parameters 



t ISM 
SPAH 
t ISM 
SMIH 

- ISM 

;w 

- ISM 

;c 



rplSM 
rplSM 



0.22 
0.11 
0.07 
0.60 
45 K 
18 K. 



(17) 



For reference, we include in Fig. [T] additional observa- 
tions of the mean Galactic cirrus emission obtained with 
Spitzer/IRAC at w avelengths 3.6, 4.5, 5.8 and 8.0 ^m 
jFlagev et al.ll2006T ) and with the COBS/Far-Infrared Ab- 
solute Spectropho tometer (FIRAS) at wavelengths between 
140 and 1000 fj,m (IDwek et al.lll997t ). The model reproduces 
these observations remarkably well, even though they were 
not included in the fit. 

We use the constraints of Fig. [T]to fix the mid-infrared 
spectral energy distribution of the ambient ISM and reduce 
the number of adjustable parameters in our model. Specif- 
ically, we fix the relative contributions by PAHs, the hot 



mid-infrared continuum and warm grains in thermal equilib- 
rium to the total infrared luminosity of the ambient ISM (i.e. 
the relative ratios of £pah> Cmir an d Cw M ) to their values of 
equation (|17|) . and we keep £cf M an d Tq SM as adjustable pa- 
rameters. Thus, for a given contribution £<i SM by cold grains 
in thermal equilibrium to the total infrared luminosity of 
the ambient ISM, the contributions by PAHs, the hot mid- 
infrared continuum and warm grains in thermal equilibrium 

g£ = 0.275(1 - £c SM ) and 
respectively. 

The total spectral energy distribution of a galaxy in our 
model is computed as the sum 

T tot _ T BC | T ISM 
^A,d — ^X,d + ^A,d i 



are £^ = 0.550(1 



^ M = 0.175(1 



c ISMn 

sc ; 



t ISM\ 



(18) 



where L^d an d are given by equations (|13[) and (|15[) . 

For some purposes, it is also convenient to define the global 
contribution by a specific dust component, including stel- 
lar birth clouds and the ambient ISM, to the total infrared 
luminosity of a galaxy. This can be written 



£ tot 
?PAH 



Cpah (1 



f lt ) + 0.550(1 



C ISM\ r 

?c ) 



(19) 
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£mir = ?mi C r (1 - U) + 0.275 (1 - C«i SM ) /„ , 
£w = £vf (1 - /„) + 0.175 (1 - £ C SM ) /„ , 

/-tot /-ISM t 



(20) 
(21) 
(22) 



for PAHs, the hot mid-infrared continuum and warm and 
cold dust in thermal equilibrium, respectively. 

We compute luminosity densities in any filter from Lx\ 
using the general formula 



C*. 



J dv v 1 L v R v 
J dvv^C v R v 



}% c J dX L\ X R\ 
c X °JdXC x \Rx 



where 



A 



/ d\ A R\ 
J dXR x 



(23) 



(24) 



is the effective wavelength of the filter of response Rx {Rv), 
and the calibration spectrum Cx (C v ) depends on the pho- 
tometric system. For t he AB system , C v is fixed at 3631 Jy, 
implying Cx oc A -2 (|Oke fc Gunnl Il983l ). For the IRAS, 
ZSO/ISOCAM and Spitzer/IRAC photometric systems, the 
convention is to use vC v = XCx = constant, and hence, 
Cx o c A" 1 (see iBeichman et all 1 19881 ; iBlommaert et all 
l2003UReach et all 120051 . and also the IRAC Data Hand- 
boolo). The Spitzer/MlPS system was calibrated using a 
blackbody spectrum of temperature 10,000 K, such that 
C x = Ba(10,000K) (MIPS Data Handbook In Fig. H 
we show an example of infrared spectral energy distribution 
L A ° d computed using our model, for parameter values typi- 
cal of normal star-forming galaxies (equation 1251 below) . As 
we shall see in Section 12.2,21 this parameterisation of in- 
frared galaxy spectra allows us to account for the full range 
of observed colours of star-forming galaxies. 



2.2.2 Comparison with observed infrared colours 

To test how well our model can reproduce the observed in- 
frared colours of galaxies in a wide range of star forma- 
tion histories, w e appeal to a samp le of 157 nearby galax- 
ies compiled by lElbaz et all (|2002f ). for whic h IRAS and 
iffO/I S OCAM observa ti ons a re available fromlrBoselli et al.l 
il 19981) . iRoussel et all (|200ll ), lLaurent et al.l i|2000l ) and 
iDale et all (|2000l ). The galaxies in this sample span wide 
ranges of morphologies, absolute infrared luminosities, 



infrared-to-blue l uminosity ratios and infrared colours (see 
lElbaz et al.ll2002l for more detail). For nearby galaxies, the 
flux density collected by the ISOCAM LW2 (6.75 (im) filter, 
-F^' 75 , tends to be dominated by the PAH emission features 
at 6.2, 7.7 and 8.6 fim. The IRAS 12-/im flux density F„ 12 
captures the PAH emission at 11.3 and 12.7 fim and the 
mid-infrared continuum emission from hot dust. The ISO- 
CAM LW3 (15 ^m) and the IRAS 25-itm flux densities, 
F^ 5 and F^ 5 , reflect primarily the mid-infrared continuum 
emission from hot dust. At longer wavelengths, the IRAS 
60- /im flux density F®° samples the emission from warm 
grains in thermal equilibrium in star-forming clouds, and 
the IRAS 100-/im flux density F* 00 that from colder grains 
in thermal equilibrium in the ambient interstellar medium. 



2 http:/ /ssc. spitzer.caltech.edu/irac/dh/ 

3 http:/ /ssc. spitzer.caltech.edu/mips/dh/ 



To test our model, we require photometric observations in 
all these filters. This reduces the sample to 107 galaxies. For 
the purpose of comparisons between our (angle-averaged) 
model and observations, we neglect possible anisotropics and 
equate flux ratios at all wavelengths to the corresponding lu- 
minosity ratios. 

In Figs. [3] and 2] we show the locations of these 
galaxies in various infrared colour-colour diagrams (black 
crosses). Fig. |3] shows Lf/Lf° as a function of jhf 
(IRAS colours), while Fig. U shows Ll 5 /L®° as a function 
of L®' 75 /L„ 5 (ISO colours). The well-known correlations be- 
tween the different infrared colours of galaxies illustrated 
by Figs. |3] and 2] suggest that the different contributors 
to the total infrared emission are related to one another. 
These relations are generally interpreted as sequences in 
the overall star formation a ctivity and dust heating (e.g. 
Hclou 19861; IDale et al.ll200ll) . Quiescent star-forming galax- 
ies with strong PAH emission and cool ambient-ISM dust 
tend to have high L^/Lf, Z^' 75 / ' h)f and jhf and 
low L®°/Lv. In contrast, actively star-forming galaxies, 
in which the mid-infrared emission is dominated by con- 
tinuum radiation by hot dust and the far-infrared emission 
by warm dust in star-forming regions, tend to have low 
Ll 2 /Lf, Ll-^/Ll 5 and hflhf and high itf/Ll 00 . 

We now use these observations to explore the influ- 
ence of each parameter of our model on the various infrared 
colours. We explore the effect of varying a single parameter 
at a time, keeping all the other parameters fixed at 'stan- 
dard' values. After some experimentation, we adopted the 
following standard parameters (corresponding to the model 
shown in Fig. [2}: 



u = 



C BC 

SPAH 

t BC 

Cmir 

£ BC 

sw 

C ISM 



rpl 



iBC 
l W 

rpISM 



0.60 
0.05 
0.15 
0.80 
0.80 
48 K 
22K. 



(25) 



These values allow the standard model to match roughly 
the observed typical (i.e. median) infrared colours of nearby 
star- forming galaxies in Figs. [3] and [4] (green circle). Each 
panel of Figs. [3] and [4] shows the effect of varying one pa- 
rameter of the model, with the other parameters held fixed 
at their standard valuesQ We can summarise the role of each 
parameter as follows. 

Fraction of total infrared luminosity contributed by dust in 
the ambient ISM. The dominant effect of increasing / M is to 
increase the contribution to the total infrared luminosity by 
cold dust. This raises the infrared luminosity around 100 /im 
(Fig. EJ, causing Lf/L™ and even Lf / 'hf to decrease 
(Figs. [3^. and|4^,). In addition, a larger / M leads to an increase 
in L^, 2 /L 25 and Zv' 75 / L^ 5 , due to the increased contribution 



4 In practice, an increase in any of £pah> £mir *>w 
accompanied by a drop in the other two fractions by virtue of 
equation 114> . Likewise, an increase in £,i SM implies a drop in the 
sum £pp^ + £)JfIR + ?w M (equation 02} . When exploring such 
variations, we keep the relative ratios of the unexplored fractions 
fixed. 



and £® c will be 
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Figure 2. Example of infrared spectral energy distribution computed using the model presented in Section 12.2,11 for 'standard' values of 
the parameters (equation [25): / p = 0.6, f^gj = 0.05, (,^ R = 0.15, = 0.80, T^ c = 48K, £,* SM = 0.8 and T£ SM = 22 K. The black 
curve shows the total infrared spectrum. The green solid curve shows the contribution by dust in the stellar birth clouds, and the green 
dashed curves the breakdown of this contribution in different components: PAHs, mid-infrared continuum and warm dust in thermal 
equilibrium. The red curve shows the contribution by dust in the ambient ISM. The total luminosity is L^ ot = 10 10 Lq (equations 1131 

ITS! and ran. 



by PAH emission features (which dominate the mid-infrared 
emission of the ISM) to the total mid-infrared emission. 

Contribution by cold dust in thermal equilibrium to the 
infrared emission from the ambient ISM. Increasing £cf M 
causes Z/^/L 25 and Z^' 75 /Lj, 5 to decrease in Figs. 03 and 
[4p, because of the corresponding drop in PAH emission from 
the ambient ISM (see footnote QJ . Also, both Ll 5 /Ll° and 
Lt° /Ll°° decrease because of the larger contribution by cold 
dust to the 60 /im flux, and even more so to the 100 /im flux. 

Contribution by PAHs to the infrared emission from stel- 
lar birth clouds. Increasing £f5St raises the contribution by 
PAH features to L®' 75 and Lj, 2 , and to a lesser extent Lj, 5 
(Fig. 0. This leads to a marked increase in L®' 75 /Lj, 5 and 
Ll 2 /Lf and a milder one in L^/Lf in Figs.EJ: andBJi. The 
slight decrease in L„° /Ll°° when £pah increases in Fig. [3]: 
is caused by the associated drop in £bq (footnote [J). 

Contribution by the hot mid-infrared continuum to the in- 
frared emission from stellar birth clouds. The main effect of 
increasing £mir i s to make I/* 5 and I/ 2 , 5 larger (Fig. [3]). This 
causes a drop in / Lf and L^/Lf, while Ll 5 /Lf in- 
creases (Figs. |3jl and |li). The slight decrease in Ll°/Ll°° 
when £mir increases in Fig. [3£ is again caused by the asso- 
ciated drop in £bq (footnote [4]). 

Contribution by warm dust in thermal equilibrium to the 
infrared emission from stellar birth clouds. Increasing ^ c 
makes Lff larger, and to a lesser extent also L 25 and L* 00 
(Fig. [2}. Remarkably, this causes L^/L 25 to drop and 



Lt°/Ll°° to rise in Fig. [3^, almost along the observational 
relation. In Fig. an increase in £^| c makes I/* 5 smaller 
(because of the associated drop in £mir! see footnote and 
hence, Ll 5 /Ll° smaller and L^ 75 /Li 5 larger. 

Equilibrium temperature of warm dust in stellar birth 
clouds. Increasing T"w C across the range from 30 to 60 K 
moves the wavelength of peak luminosity of dust in thermal 
equilibrium in the stellar birth clouds roughly from 70 to 
40 fim. This causes Li 2 /L 25 to drop and Ll°/Ll°° to rise in 
Fig. if. In Fig.gf, the effects on L^ 75 / Ll 5 and Ll 5 /Ll° are 
negligible. 

Equilibrium temperature of cold dust in the ambient ISM 
(not shown). The most significant effect of increasing Tq SM 
at fixed £c SM in Figs. [3] and g] is a small rise in L 6 °/Lf° 
caused by a blue shift of the peak infrared luminosity of cold 
dust in the ambient ISM (from about 140 to 80 (im as Tq SM 
increases from 15 to 25 K). Variations in Tq SM also have a 
significant influence on the emission redward of 100 /im. 

The green lines in Figs. [3] and g] show that the extrem- 
ities of the observational relations between the different in- 
frared colours of galaxies cannot be reached by varying a 
single model parameter at a time. This suggests that varia- 
tions in the different dust components of galaxies are related 
to each other. However, we have checked that the properties 
of every galaxy in Figs. [3] and g] could be reproduced with 
at least one combination of parameters of our model. We 
illustrate this by showing two models lying at the ends of 
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Figure 3. Ratio of 60-fim to 100-/^m IRAS lumino sity density p lotted against ratio of 12-/^m to 25-fim IRAS luminosity density. The 
data points in all panels (black crosses) are from the lElbaz et al.l J2002T) sample discussed in Section 12 21 (typical measurement errors are 
indicated in the upper left panel). In each panel, the green line shows the effect of varying one parameter of the model from the lower 
end of the range (square) to the standard value (green circle) to the upper end of the range (triangle), with all other parameters fixed 
at their standard values: (a) fraction of total infrared luminosity contributed by dust in the ambient ISM, = 0.05, 0.50 and 0.95; (b) 
contribution by cold dust in thermal equilibrium to the infrared luminosity of the ambient ISM, f;,i SM = 0.50, 0.80, 1.0; (c) contribution by 
PAHs to the infrared luminosity of stellar birth clouds, = 0.00, 0.05 and 0.50; (d) contribution by the hot mid-infrared continuum 

to the infrared luminosity of stellar birth clouds, £mir = 0.00, 0.15 and 0.50; (e) contribution by warm dust in thermal equilibrium to 
the infrared luminosity of stellar birth clouds, £w = 0.15, 0.80 and 0.95; (f) equilibrium temperature of warm dust in stellar birth 
clouds, T^ c = 30, 48 and 60 K. In all panels, the red circle corresponds to the same quiescent model galaxy dominated by cold dust 
emission, while the blue circle corresponds to an actively star-forming galaxy dominated hot dust emission (sec Table [T] and Section l2.2l 
for a description of these models). 



the observational relations; a 'cold' infrared spectrum char- 
acteristic of a quiescent galaxy with little star formation 
(red circle); and a 'hot' infrared spectrum characteristic of 
an actively star- forming, starburst galaxy (blue circle). The 
parameters of these models are listed in Table [1] We empha- 
sise that these are not unique sets of parameters optimised 
to fit specific galaxy spectral energy distributions. Rather, 
they are examples of how the colours in those regions of the 
diagrams can be reproduced using our model. 

Figs. [5] and H] allow us to draw some general conclusions 
about the influence of the various parameters of our model 
on the observed infrared colours of galaxies. For example, 
as expected from Fig. [2j the L®°/Ll°° colour appears to be 
controlled primarily by the fraction / M of total infrared lu- 
minosity contributed by dust in the ambient ISM and the 
properties of dust in thermal equilibrium: the relative con- 



tribution £ c by cold dust to the infrared luminosity of the 
ambient ISM (and the temperature Tq SM of this dust ) and 
the relative contribution by warm dust to the infrared 
luminosity of stellar birth clouds (and the temperature T^ c 
of this dust). These parameters also have distinct effects on 
the mid-infrared colours Lj 2 / 'Lf , Ll^/Lf and L^/Lf, 
indicating that they can be constrained independently from 
fits of extended infrared spectral energy distributions. The 
mid-infrared colours are primarily controlled by the different 
components of hot dust: the relative contributions £f^h an d 
£mir °f PAHs and the hot mid-infrared continuum to the 
the infrared luminosity of stellar birth clouds, and the con- 
tribution by PAHs to the infrared luminosity of the ambient 
ISM, which is controlled indirectly by £c SM . I n Section[3]be- 
low, we show how well these various model parameters can 
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Figure 4. Ratio of 15-fim to 60-/^m /SO and IRAS luminosity densit y plott ed against ratio of 6.75-/xm to 15-/im ISO luminosity density. 
The data points in all panels (black crosses) are from the lEIbaz et alj l l2002h sample discussed in Section 12.21 In each panel, the green line 
shows the effect of varying one parameter of the model from the lower end of the range (square) to the standard value (green circle) to 
the upper end of the range (triangle), with all other parameters fixed at their standard values. The models are the same as in Fig. [3] In all 
panels, the red circle corresponds to the same quiescent model galaxy dominated by cold dust emission, while the blue circle corresponds 
to an actively star-forming galaxy dominated hot dust emission (see Table [T] and Section l2.2l for a description of these models). 



Table 1. Parameters of the 'standard', 'cold' and 'hot' models 
plotted in Figs. [3] and [4] 



Model 


'cold' 


'standard' 


'hot' 


U 


0.75 


0.60 


0.20 


fISM 


0.75 


0.80 


0.90 


cBC 
?PAH 

cBC 
?MIR 

;BC 

?w 


0.45 
0.15 
0.40 


0.05 
0.15 
0.80 


0.01 
0.09 
0.95 


T* C (K) 
T(J SM (K) 


10 
18 


18 
22 


55 
25 



2.2.3 Constraints on dust mass 

It is also of interest to derive constraints on the dust mass in 
galaxies. The mass Md(Td) in dust grains in thermal equi- 
librium at the temperature Td can be estimated from the 
far-infrared radiat ion L^ d of these grains using the formula 
l|Hildebrandlll983l ) 

Ll d = 4tt M d (T d ) k\ B\(Td) , (26) 

where n\ and B\ have been defined before (equation II Op . 
We adopt this formula to estimate the mass contributed 
by dust in thermal equilibrium in stellar birth clouds (with 
temperature Td = 7w°) an d in the ambient ISM (with tem- 



peratures T d = T^ SM 



and T d = T<f M ) in our model. We 



adopt as before a dust emissivity index /3 = 1.5 for warm 
dust and (3 = 2 for cold dust an d normalise k,\ at 850 fim 
assuming kssom™ — 0.77g _1 cm 2 l|Dunne et alj|200(ih . Using 
equations l|10|l. (fT3|) . (fT5|) and ([26l . we compute the mass 
contributed by warm dust in stellar birth clouds and in the 
ambient ISM as 



be constrained in galaxies with observed infrared spectral 
energy distributions. 



Sw 1 



l-/„)£d 



4tt 



cI\k x B X (T* C ) 



(27) 
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and 



M4 



C ISM r 

tw Jm 



4tt 



(28) 



and that contributed by cold dust in the ambient ISM as 



]\yf ISM C ISM r 

M c — ?c 7m 



4tt 



(29) 



To include the contribution by stochastically heated 
dust grains (not in th ermal equilibr i um; se e Section [2.2.1[) , 
we adopt a standard iMathis et al.l (1 19771 ) distribution of 
grain sizes N(a) oc a -3 ' 5 over the range 0.005 fim < a < 
0.25 /im. We assume that the stochastically heated dust 
grains are very small (a < 0.01 fim) and have mass densities 
typical of graphite, p ~ 2.26gcm , while bigger grains in 
thermal equilibriu m have mass densiti es typical of silicates, 
p ~ 3.30 gem" " 3 ijDraine fc Leelll984r i. For these assump- 
tions, the mass contributed by grains of all sizes is about 
5 per cent larger than that contributed by big grains alone 
(a > 0.01 pm). The contribution by PAHs to the overall 
dust mass is also ve ry small, of the order of a few per cent 
l|Draine et alj|2007h . We therefore estimate the total dust 
mass of a galaxy as 



M A » 1.1 (M^ c + M^ SM + M£ s 



(30) 



where M^ c , and M(i SM are given by equations (J27J), 

(l28l) and (l29l. 



2.3 Combined ultraviolet, optical and infrared 
spectral energy distributions 

A main feature of our model is the consistent modelling 
of ultraviolet, optical and infrared spectral energy distribu- 
tions of galaxies. This is achieved by first computing the 
total energy absorbed by dust in stellar birth clouds and in 
the ambient ISM, hf° and hf M (Section EU, and then 
re-distributing it at infrared wavelengths (Section 12. 2[) . The 
main assumptions are the conservation of the energy ab- 
sorbed and reradiated by dust, and that the dust in the 
ISM of galaxies is heated only by starlight (in particular, we 
ignore the possible influence of an active galactic nucleus). 
Different combinations of star formation histories, metal- 
licities and dust contents can lead to the same absorbed 
energies L^ 30 and L^ SM in a model galaxy. Furthermore, 
these energies can be distributed in wavelength using dif- 
ferent combinations of dust parameters in the stellar birth 
clouds (£pahi £mir> Cw° an d Jw°) and the ambient ISM 
(C C SM and T^ SM ). In our model, therefore, a wide range of 
ultraviolet and optical spectral energy distributions can be 
associated to a wide range of infrared spectral energy dis- 
tributions, at fixed L^ c and Lj SM (or equivalently, at fixed 
/ M and L d ot ; see equations [5] and [8} . In Section [3] below, we 
show how combined observations at ultraviolet, optical and 
infrared wavelengths can be used to uniquely constrain the 
star formation histories and dust properties of galaxies using 
this model. 

To illustrate the combination of ultraviolet, optical and 
infrared spectral energy distributions with our model, we 
compute examples of attenuated stellar population spectra 
consistent with the / M parameters of the cold, standard and 
hot infrared models of Table[T1(Section l2.2.2|) . For simplicity, 



we assume solar metallicity and select models with exponen- 
tially declining star formation rates 



i>{t) oc exp(-jt) . 



(31) 



where 7 is the star formation timescale parameter. We 
choose models with 7 = 0, 0.07 and 0.25 Gyr -1 at the ages 
t = 1.4, 10 and 10 Gyr, respectively, to represent a starburst, 
a no rmal star-forming and a quie scent star-forming galaxy 
(e.g. iKauffmann fc Charlotl [l998h . For the attenuation of 
starlight by dust (equations [3] and |4}, we adopt effective 
dust absorption optical depths Tv = 2.0, 1.5 and 1.0, for the 
starburst, normal star-forming and quiescent star-forming 
models, respectively. These values are consistent with the ex- 
pectation that more actively star- form i ng galaxies are more 
obscured (e.g. Wang fc Heckmanll 19961 ; lllopkins et al. I l200ll ; 

200 11 ). We adopt fractions p — 0.1, 0.3 and 



ISullivan et al.l 



0.5 of fy arising from dust in the ambient ISM. The result- 
ing model galaxies have / M ~ 0.2, 0.6 and 0.75, respectively, 
consistent with the parameters of the hot, standard and cold 
infrared models of Table [T] 

In Fig. [S] we show the combined spectral energy distri- 
butions of the starburst + hot infrared, normal star-forming 
+ standard infrared, and quiescent star-forming + cold in- 
frared models, after scaling in each case the infrared lumi- 
nosity to the total luminosity L^ ot absorbed by dust. As we 
shall see in the next section, the ability to compute such 
combined ultraviolet, optical and infrared spectral energy 
distributions for wide ranges of physical parameters of galax- 
ies has important implications for statistical estimates of 
star formation histories and dust properties. 



3 CONSTRAINTS ON PHYSICAL 
PARAMETERS FROM 
MULTI- WAVELENGTH GALAXY 
OBSERVATIONS 

In this section, we show how the simple model described 
in Section [2] can be used to extract star formation histories 
and dust properties from ultraviolet, optical and infrared 
observations of galaxies. In Section T3. II we first describe our 
methodology to derive statistical constraints on galaxy phys- 
ical parameters from multi-wavelength observations. Then, 
in Section 13.21 we use this methodology to constrain the 
physical parameters of star-forming galaxies in the SINGS 
sample. We compare our results to those that would be ob- 
tained using previous models in Section f3. 31 Finally, we dis- 
cuss possible sources of systematic errors associated with our 
approach in Section f3. 41 and summarize the applicability of 
this model in Section 1331 



3.1 Methodology 

The model of Section [5] allows one to compute the ultravio- 
let, optical and infrared emission from galaxies. This model 
contains a minimum number of adjustable parameters re- 
quired to account for the observed relations between various 
integrated spectral properties of galaxies: age, star formation 
history, stellar metallicity, two components for the attenua- 
tion by dust and four contributors to the infrared emission 
(PAHs, hot mid-infrared continuum, warm and cold dust 
in thermal equilibrium). The large number of observable 
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Figure 5. Examples of spectral energy distributions obtained by combining the infrared models of Table [T] with attenuated stellar 
population spectra corresponding to the same contributions by dust in stellar birth clouds (1 — f^) and in the ambient ISM (/^) to the 
total energy L^ ot absorbed and reradiated by dust (Section 12.31 . (a) Quiescent star-forming galaxy spectrum combined with the 'cold' 
infrared model of TablelTI (b) normal star-forming galaxy spectrum combined with the 'standard' infrared model of Tablefl] (c) starburst 
galaxy spectrum combined with the 'hot' infrared model of Table IT1 (see text for details about the parameters of the stellar population 
models). Each panel shows the unattenuated stellar spectrum (blue line), the emission by dust in stellar birth clouds (green line), the 
emission by dust in the ambient ISM (red line) and the total emission from the galaxy, corresponding to the sum of the attenuated stellar 
spectrum and the total infrared emission (black line). 



quantities to which the model can be compared insures that 
these different adjustable parameters can be constrained in 
a meaningful way (see below). A usual limitation of this type 
of study is that several different combinations of physical pa- 
rameters can lead to similar spectral energy distributions of 
galaxies. For example, age, metallicity and dust attenuation 
have similar effects on the ultraviolet and optical spectra of 
galaxies. An efficient way to derive statistical constraints on 
the various parameters in these conditions is to consider a 
wide library of models encompassing all plausible parameter 
combinations. Given an observed galaxy, we can build the 
likelihood distribution of any physical parameter by evalu- 
ating how well each model in the library can account for the 
observed properties of the galaxy. This B ayesian approach 
is similar to that used, for example, by iKauffmann et abl 



(2003) to interpret the optical spectra of SDSS galaxies. 
The underlying assumption is that the library of models is 
the distribution from which the data were randomly drawn. 
Thus, the prior distribution of models must be such that 
the entire observational space is reasonably well sampled, 
and that no a priori implausible corner of parameter space 
accounts for a large fraction of the models. 

3.1.1 Model library 

We build a comprehensive library of models by generating 
separately a random library of stellar population models, 
for wide ranges of star formation histories, metallicities and 
dust contents, and a random library of infrared spectra, for 
wide ranges of dust temperatures and fractional contribu- 
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tions by the different dust components to the totai infrared 
luminosity. We then combine these libraries following the 
procedure outline in Section [2.31 

For simplicity, we follow iKauffmann et all |2003) and 
parametrise each star formation history in the stellar pop- 
ulation library in terms of two components: an underlying 
continuous model, characterised by an age t s and a star for- 
mation time-scale parameter 7 (equation I31f) . and random 
bursts superimposed on this continuous model. We take t s 
to be uniformly distributed over the interval from 0.1 and 
13.5 Gyr. To avoid oversampling galaxies with negligible 
current star formation, we distribute 7 using the probabil- 
ity density function p(~y) = 1 — tanh(87 — 6), which is ap- 
proximately uniform over the interval from to 0.6 Gyr -1 
and drops exponentially to zero around 7 = 1 Gyr -1 . Ran- 
dom bursts occur with equal probability at all times un- 
til i g . We set the probability so that 50 per cent of the 
galaxies in the library have experienced a burst in the 
past 2 Gyr. We parametrise the amplitude of each burst 
as A = Mburst/A^cont, where Mb mst is the mass of stars 
formed in the burst and M cont is the total mass of stars 
formed by the continuous model over the time t s . This ratio 
is distributed logarithmically between 0.03 and 4.0. During 
a burst, stars form at a constant rate over the time tburst, 
which we distribute uniformly between 3 x 10 7 and 3 x 10 s yr. 
We distribute the models uniformly in metallicity between 
0.02 and 2 times solar. 

We sample attenuation by dust in the library by ran- 
domly drawing the total effective V-band absorption opti- 
cal depth, fv, and the fraction of this contributed by dust 
in the ambient ISM, fi (equations [3] and 3J). We distribute 
tv according to the probability density function p(fv) = 
1 — tanh(1.5fv — 6.7), which is approximately uniform over 
the interval from to 4 and drops exponentially to zero 
around tv = 6. For fi, we adopt the same probability density 
function as for 7 above, i.e. p(p) ~ 1 — tanh(8 fi — 6). We note 
that these priors for attenuation encompass the dust prop- 
erties of SDSS galaxies, for which fv and 11 peak around 1.0 
and 0.3, respectively, w ith broad scatter |Brinchmann et all 
l2004lKong et al.l 120041 ). Our final stellar population library 
consists of 50,000 different models. 

In parallel, we generate a random library of infrared 
spectra as follows. We take the fraction / M of the total in- 
frared luminosity contributed by dust in the ambient ISM 
to be uniformly distributed over the interval from to 1. We 
adopt a similar distribution for the fractional contribution 
by warm dust in thermal equilibrium to the infrared lumi- 
nosity of stellar birth clouds, For each random drawing 
of £w°j we successively draw the contributions by the other 
dust components to the infrared luminosity of stellar birth 
clouds (i.e., hot mid-infrared continuum and PAHs) to sat- 
isfy the condition in equation (| 14f) : we draw £mir from a 
uniform distribution between and 1 — £, BC , and we set 



BC 



£mir- While this procedure does not ex- 



clude values of £mir an d Crah close to unity, it does favour 
small values of these parameters, and hence, it avoids over- 
sampling physically implausible models. We take the equi- 
librium temperature T^ c of warm dust in the stellar birth 
clouds to be uniformly distributed between 30 and 60 K, and 
that Tq SM of cold dust in the ambient ISM to be unifomly 
distributed between 15 and 25 K. We draw the fractional 
contribution £cf M by cold dust in thermal equilibrium to 



the infrared luminosity of the ambient ISM from a uniform 
distribution between 0.5 and 1 (this also defines the contri- 
butions £p!h, Cmir and by PAHs, the hot mid-infrared 
continuum and warm dust to the infrared luminosity of the 
ambient ISM, as described in Section ^, 2. l[l . Our final library 
of infrared spectra consists of 50,000 different models. 

We combine the library of stellar population models and 
that of infrared spectra by associating together models with 
similar / M , which we scale according to the total infrared 
luminosity L^ ot , as outlined in Section 12.31 In practice, we 
associate each model in the stellar population library to all 
the models in the infrared spectral library that have similar 
fn to within some error interval 5/^0 That is, we associate 
each stellar population spectrum, characterised by / M , 
to all the infrared spectra characterised by f/j, m , such that 
/ M = /^ FH i Sfp. Each spectral combination satisfying 
this condition is included in the final model library and is as- 
signed a value = (/ M SFH + / M IR )/2. We adopt 5/ M = 0.15, 
which allows good reproductions of combined ultraviolet, 
optical and infrared observations of galaxies ( Section 13. 2[) . 

Our final library of combined ultraviolet, optical and 
infrared spectral energy distributions consists of about 661 
million models. It is important to note that such a large 
number of models is required to properly sample the multi- 
dimensional observational space. 



3.1.2 Statistical constraints on physical parameters 

We now investigate the accuracy to which the model library 
described in the previous section can help us constrain the 
star formation histories and dust properties of galaxies for 
which multi- wavelength observations are available. To assess 
this, we evaluate how well we can recover the parameters of a 
random set of models with known properties, based on spec- 
tral fits. We consider the following observable quantities: the 
GALEX far-ultraviolet (FUV, 1520 A) and near-ultraviolet 
(NUV, 2310 A) luminosities; the optical UB V luminosities; 
the 2MASS near-infrared JHKs (1.25, 1.65 and 2.17 (j,m) lu- 
minosities; the Spiteer/IRAC 3.6, 4.5, 5.8 and 8.0 fim lumi- 
nosities; the ISO/ISOCAM 6.75 and 15 /im luminosities; the 
IRAS 12, 25, 60 an d 100 /im luminosities; the Spitzer/MLPS 
jRieke et al.ll2004l) 24, 70 and 160 luminosities; and the 
SCUBA jHolland et all 119991 ) 850 luminosity. We also 
include the hydrogen Ha and H/3 recombination-line lumi- 
nosities. 

We compute these observable quantities from the spec- 
tra of all 661 million models in the library and randomly 
select a subset of 100 models. To mimic observational con- 
ditions, we perturb the luminosities of this subset of models 
assuming a fixed uncertainty of 10 per cent. We then per- 
form spectral fits to recover the likelihood distributions of 
the physical parameters of these 'mock galaxies' as follows. 
We first compare the luminosities of each mock galaxy in 
the sample to the luminosities of every model j in the li- 
brary to measure the v 2 goo dness-of-fit of that model (e.g., 



iBevington fc Robinsonj|2003l) 



5 This allows us to include the uncertainties that could arise, for 
example, from orientation effects, in the connection between the 
stellar and dust emission. 
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Figure 6. Median-likelihood estimates of 12 parameters (indicated in the upper left corner of each panel) recovered from the spectral fits 
of 100 mock galaxies, compared to the true values of these parameters (see Section 13.1.20 . The error bars indicate the 16-84 percentile 
ranges in the recovered probability distributions of the parameters. These results were obtained by fitting simulated GALEX (FUV and 
NUV), optical ( UBV), 2MASS (JHKs), Spitzer/IRAC (3.6, 4.5, 5.8 and 8.0 ftm), KO/ISOCAM (6.75 and 15ftm), IRAS (12, 25, 60, 
100 lira), Spitzer/MlPS (24, 70, and 160 (im). SCUBA (850 ftm), and Ha and H/3 luminosities. 



2 



(32) 



where L l v and L„ j are the luminosities in the i band of 
the mock galaxy and the j th model, respectively, en is the 
(10 percent) uncertainty in L\,, and 



T % 



(33) 



is the model scaling factor that minimizes Xj - Then, we build 
the probability density function of any physical parameter 
of the mock galaxy by weighting the value of that param- 
eter in the j th model by the probability exp(— Xj/2). We 
take our 'best estimate' of the parameter to be the median 
of the resulting probability density function and the asso- 
ciated confidence interval to be the 16-84 percentile range 



(this would be equivalent to the ± la range in the case of a 
Gaussian distribution). 

In Fig. [6] we compare the median-likelihood estimates 
to the true values of 12 parameters recovered in this way 
for the 100 mock galaxies in our sample. These parameters 
are: the fraction of total infrared luminosity contributed by 
dust in the ambient ISM, / M ; the total effective V-band ab- 
sorption optical depth of the dust, tv> the fraction of this 
contributed by dust in the ambient ISM, the specific star 
formation rate, i/>g, defined as the ratio 



Mt) = 



(34) 



t 8 M, 

of the star formation rate averaged over the past t& = 10 8 yr 
to the current stellar mass M* of the galaxy; the stellar 
mass, M* (this accounts for the mass returned to the ISM 
by evolved stars); the total infrared luminosity, L^ Qt ; the 
fractional contributions by PAHs, the hot mid-infrared con- 
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tinuum and warm dust in thermal equilibrium to the in- 
frared luminosity of stellar birth clouds, £pahj £mir an d 
£^ c ; the equilibrium temperature of warm grains in stel- 
lar birth clouds, T^°; the contribution by cold dust to the 
total infrared luminosity of the ambient ISM, £cf M , an d the 
equilibrium temperature of this dust, Tq Sm . Most of these 
parameters are recovered remarkably well by our model. In 
particular, / M , f v , M , </>s, M„ if, £mir, & sm and T^ SM 
are recovered to very high accuracy. The most uncertain 
parameters are £pah> £w° and T^ c , for which the typical 
confidence intervals reach almost 0.14, 0.18 and 11 K, respec- 
tively. We conclude that our model is a valuable tool for de- 
riving statistical constraints on the star formation histories 
and dust properties of galaxies for which multi-wavelength 
(ultraviolet, optical and infrared) observations are available. 
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3.2 Application to an observational sample 

3.2.1 The sample 

Here, we exploit our model to interpret a wide range of ul- 
traviolet, optical and infrared observations of a sample of 
well-studied nearby gal axies: the Spitzer Infr ared Nearby 
Galaxy Survey (SINGS; iKennicutt et alj|2003h . This sam- 
ple contains 75 galaxies at a median distance of 9.5 Mpc 
(we adopt a Hubble constant Ho = 70kms _1 Mpc _1 ). The 
galaxies span wide ranges in morphology (from E-S0 to Im- 
10) and star formation activity (from quiescent to starburst). 
Some galaxies include low-luminosity active galactic nuclei 
(AGNs)0 We note that this sam ple extends out to lower 
total infrared luminosities than the lElbaz et al.l (|2002T > sam- 
ple used to calibrate the infrared properties of our model 
in Section (L d tot < lO 11 !,© instead of L\ ot < 1O 12 L ), 

but it includes observations across a much wider range of 
wavelengths. 

Observations at ultraviolet, optical and infrared wave- 
lengths are available for most galaxies in this sample. 
GALEX ultraviolet (NU V and FUV) ob servations are 
available for 70 galaxies l|Dale et al.l l2007f l. In the opti- 
cal, we adopt UBV fluxes from the RC3 for 65 galaxies 
|de Vaucouleurs et a l. 1991 )Q Near- infrared JHKs fluxes are 
avai lable for all galaxies from the 2MASS Large Galaxy At- 
las (jjarrett et al.l 120031 ) . In the mid- a nd far-infra r ed, w e 
use the Spitzer observations published bv lDale et aill|2007h . 
This includes IRAC photometry at 3.6, 4.5, 5.8 and 8.0 jim. 
and MIPS photometry at 24, 70 and 160 pm. ISO, IRAS 
and SCUBA 850 /im observations are also available for 13, 
65 and 22 SINGS galaxies, respectively. These observations 
are typically of lower quality than Spitzer data, with photo- 
metric uncertainties around 20 per cent for ISO and IRAS 
and 30 per cent for SCUBA. More detail about the photom- 
etry of SINGS galaxies can be found in Dale et al. I (|2005l . 
l2007l ). 



6 The low-luminosity AGNs should have a negligible impact on 
the integrated broad-band fluxes of the galaxies. We distinguish 
galaxies hosting low-luminosity AGNs from normal star-forming 
galaxies in the spectral analyses presented later in this paper. 

7 We choose not to use the optical BVRI fluxes tabulated by 
Dale et al. 1 1120071 ). because of calibration problems in the optical 
photometry of the SINGS Fifth Enhanced Data Release (D. Dale 
& A. Gil de Paz, priv. comm.) 
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Figure 7. Selected properties of the SINGS galaxy sample dis- 
cussed in Section 13.21 (a) Ratio of 70-fira to 160-/xm luminosity 
plotted against ratio of 8-fim to 24-/im luminosity, (b) Ratio of 
total-infrared to ultraviolet luminosity as a function of ultraviolet 
spectral slope (the I RX— U V diagram) . The dotted line in this dia- 
gram shows a fit byjKong ii ei ii alJj200J) to the relation followed by 
starburst galaxies llMeurer et al.1 19991 ) . In both p anels, crosses re- 
fer to galaxies with no AGN activity according to Kennicutt et al. 
(2003), while open squares refer to galaxies hosting low-luminosity 
AGNs (LINER and Seyfert). Filled circled indicate three galaxies 
selected along the L^/Lj, 60 versus Z/*/L 24 relation, for which 
a wide set of observations are available (Section I3.2.2H . Typical 
measurement errors are indicated in each diagram. 



It is important to note that these multi-wavelength 
observations can be combined in a meaningful way, since 
the fluxes all pertain to the total emission from a galaxy. 
For example, the GALEX (FUV and NUV), RC3 (UBV) 
and 2MASS (JHKs) fluxes were obtained by integrating 
extrapolated surface brightness profiles to include the en- 
tire emission from a galaxy (seelde Vaucouleurs et al.lll99ll ; 
IJarrett et all 120031 ; iGil de Paz et all 120061 ). Also, the in- 
frared Spitzer IRAC and MIPS fluxes include extended- 
source ap erture corrections (amo unting typic ally to 10 per 
cent; see iDale et al.l 120051 . l2007f ). Following iDraine et aH 
l|2007h . we exclude 9 galaxies with bad IRAC and MIPS 
detections (because of contamination by external sources; 
very low surface brightness compared to foreground cir- 
rus emission; and saturation issues): NGC 584, NGC 3034, 
NGC 1404, NGC 4552, M81 DwarfA, M81 DwarfB, 
DDO 154, DDO 165 and Holmberg IX. For 19 galaxies, we 
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include Ha and H/9 emission-line fluxes (corrected for con- 
tamination by stellar absorption) from the integrated spec- 
troscopy of iMoustakas fc Kennicuttl (|2006h . Our final sam- 
ple includes 66 galaxies, of which 61 have GALEX measure- 
ments. As before, for the purpose of comparisons with mod- 
els, we neglect possible anisotropies and equate flux ratios 
at all wavelengths to the corresponding luminosity ratios. 

We illustrate the ultraviolet and infrared properties 
of the SINGS galaxies in Fig. [7] Fig. G$, shows the MIPS 
Ll°/Ll 60 luminosity ratio plotted against the IRAC+MIPS 
Lt/L„ 4 luminosity ratio. SINGS galaxies in this diagram 
follow a se quence very s imilar to that followed by the galax- 
ies of the lElbaz et al.l (|2002l ) sample of Section 12.21 in the 
analogous IRASLf/Ll 00 versus hfjhf diagram (Fig.©. 
In Fig. [7p, we plot ratio of total-infrared to ultraviolet lu- 
minosity L*j° / L-fuv as a function of ultraviolet spectral 
slo pe ffijy (the 'IRX-U V diagram'). We used equation (4) 
of bale fc Heloul <|2002i ) to estimate L^ ot from the MIPS 
observation s at 2 4, 70 and 160 /im, and equation (1) of 
iKong et al.l l|2004h to compute f3 uv from the GALEX FUV 
and NUV luminosities. Previous studies have shown that 
starburst galaxies follow a tight sequence in this diagram 
(indicated by the dotted line in Fig. [Tl 3 ), while more qui- 
escent star-forming galaxies tend to exhibit redder ultra- 
violet spectra (i. e. larger Pvv) than starburst galaxies at 
fixed L\ ot /L FUV l|Meurer et al.lll999l ; [Belill2002l ; IKong et all 
|2004| ). Fig. [3 confirms that the SINGS galaxies are repre- 
sentative of the ultraviolet and infrared properties of nearby 
galaxies. 

3.2.2 Example of constraints on physical parameters 

We now use the methodology outlined in Section [37T]to con- 
strain the star formation histories and dust properties of the 
SINGS galaxies. We first assess how well our model can re- 
produce the observed multi- wavelength observations of these 
galaxies. For each observed galaxy, we select the model of 
the library presented in Section l3TT1 that minimizes Xj (equa- 
tion [32]), as computed by including as many luminosities as 
available in the following bands: GALEX (FUV and NUV), 
optical (UBV), 2MASS (JHKs), Spitzer/IRAC (3.6, 4.5, 5.8 
and 8.0 /im), /SO/ISOCAM (6.75 and 15/xm), IRAS (12, 
25, 60 and 100 /mi), Spitzer/MlPS (24, 70 and 160 /im), 
SCUBA (850 /im) and Ha and H^. Fig. [8] shows the re- 
sulting distribution of the difference between the observed 
luminosity L° hB and the best-fit model luminosity L™ od , in 
units of the observational error <r obs , for the 66 galaxies in 
the sample (plain histograms). Each panel shows the result 
for a different photometric band. For reference, the shaded 
histograms show the results obtained when excluding the 
22 galaxies iden t ified as hosts of low-luminosity AGNs by 
iKennicutt et al l (|2003l ). 

Fig. [8] shows that our model can reproduce simulta- 
neously the ultraviolet, optical and near-, mid- and far- 
infrared observations of all but a few SINGS galaxies to 
within the observational errors. This is remarkable, consid- 
ering the relative simplicity of the model. We find that the 
galaxies the least well fitted at ultraviolet and optical wave- 
lengths tend to be low-metallicity dwarf gal axies, such as the 
blue compact dwarf galaxy NGC 2915 (e.g. iLee et al.ll2003l ) 
and the diffuse dw arf irregular galaxy Holmberg II (e.g. 
iHoessel et al.lfl998l ). The difficulty in reproducing these ob- 



servations may arise either from a limitation of the spectral 
synthesis code, or of the simple dust prescription, or from 
an underestimate of observational errors. The three galaxies 
for which the models significantly underestimate the 3.6 /tin 
emission in Fig. [8] are also low-metallicity dwarf galax- 
ies (NGC4236, IC4710 and NGC6822). For these galaxies, 
our simple scaling of the near-infrared continuum strength 
around 4 fim with the flux density of the 7.7 /mi PAH feature 
may not be appropriate (see Section 12. 2. 1 [I . Fig. [H] further 
shows that the galaxies with 8 /mi mid-infrared luminosities 
least well reproduced by the model tend to be hosts of low- 
luminosity AGNs. In these galaxies, the contribution to the 
infrared emission by dust heated by the AGN may be sig- 
nificant. Despite these few outliers, we conclude from Fig. [8] 
that our model is adequate to investigate the constraints set 
on the physical properties of SINGS galaxies by their ultra- 
violet, optical and infrared spectral energy distributions. 

To investigate the implications of these fits for the deter- 
mination of physical parameters, we focus on three galaxies 
spanning different dust properties along the LlP/L^ 60 ver- 
sus L%jL 24 colour-colour relation in Fig. [7^., for which a 
wide set of observations are available: NGC 3521, a spiral 
(SABbc) galaxy at a distance of 9 Mpc; NGC 337, a spiral 
(SBd) galaxy at a distance of 24.7 Mpc; and Mrk 33, a dwarf 
irregular (Im) starburst galaxy at a distance of 21.7 Mpc. 
For all three galaxies, observations from GALEX (FUV and 
NUV), 2MASS (JHKs), Spitzer (IRAC and MIPS), IRAS 
and SCUBA are avail able, along with integrated Ha and 
H/3 spectroscopy from IMoustakas fc Kennicuttl (|2006l ). For 
NGC 3521 and NGC 337, additional constraints are avail- 
able from RC3 ( UBV) and ISO (ISOCAM 6.75 and 15/tm) 
photometry. 

Fig. [9] shows the models (in black) providing the best 
fits to the observed ultraviolet, optical, Ha and H/3, and in- 
frared spectral energy distributions of these three galaxies 
(in red). The quality of the fits is remarkable, as quanti- 
fied by the residuals shown at the bottom of each spectrum. 
We note that, according to the best-fit models, the sequence 
of increasing F™/F™ and decreasing F^/F 24 colours from 
NGC 3521 to NGC 3521 to Mrk 33, which reflects a drop 
in the relative intensity of PAHs and a blueshift of the peak 
infrared luminosity (i.e. a rise in the overall dust tempera- 
ture), is associated to an increase in star formation activ- 
ity. This is apparent from the relative strengthening of the 
unattenuated ultraviolet and optical spectrum (in blue) from 
NGC 3521 to NGC 3521 to Mrk 33 in Fig. [3 We investigate 
this trend further in Section 13.2.31 below. 

We can study in more detail the constraints set on the 
star formation histories and dust properties of the galax- 
ies by looking at the probability distributions of the corre- 
sponding model parameters. In Figs. [TU] and 111! we show 
the likelihood distributions of 14 quantities constrained by 
the observed spectral energy distributions of NGC 3521, 
NGC 337 and Mrk 33. These distributions were derived us- 
ing the model library of Section [3.1.11 following the proce- 
dure outlined in Section 13.1.21 Fig. [10] shows the likelihood 
distributions of parameters pertaining to the star formation 
history and dust content: fraction of the total infrared lumi- 
nosity contributed by dust in the ambient ISM (/ M ); fraction 
of the total effective l/-band absorption optical depth of the 
dust contributed by the ambient ISM (/x); effective l/-band 
absorption optical depth of the dust (fy); effective V-band 
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Figure 8. Distribution of the difference between observed luminosity L^ s and best-fit model luminosity L™ od , in units of the observa- 
tional error cr obs , for the 66 galaxies of the SINGS sample studied in Section 13,21 Each panel refers to a different photometric band (as 
indicated). The best-fit model for each galaxy was selected by fitting as many luminosities as available in the following bands: GALEX 
{FUVand NUV), optical (UBV), 2MASS (JHKs), Spitzer/IRAC (3.6, 4.5, 5.8 and 8.0 /xm), JSO/ISOCAM (6.75 and 15/mi), IRAS (12, 
25, 60 and 100 /mi), Spitzer/MTPS (24, 70 and 160 /mi), SCUBA (850 /on) and Ha and H/3. In each panel, the plain histogram shows 
the result for the full sample of 66 SINGS galaxies, while the shaded histogram shows the result for the subsample of 44 galaxies with 
no AGN activity according to Kenn icutt et alJ l|2003h . The dotted curve shows a Gaussian distribution with unit standard deviation, for 
reference. 



absorption optical depth of the dust in the ambient ISM 
(/ify); specific star formation rate {ij)s\ equation I34[); stel- 
lar mass (M t ); and dust mass (M^; equation I30[) . Fig. 1111 
shows the likelihood distributions of parameters pertaining 
to the dust emission: total infrared luminosity (L^ ot ); global 
contributions (i.e. including stellar birth clouds and the am- 
bient ISM) by PAHs (£pah)j the hot mid-infrared continuum 
(£m?r) and warm dust in thermal equilibrium to the 

total infrared luminosity (equations 1191 1201 and I21|l ; equi- 
librium temperature of warm dust in stellar birth clouds 
(Tw°); contribution by cold dust in thermal equilibrium to 
the total infrared luminosity (£rj ot > equation[22} ; and equilib- 
rium temperature of cold dust in the ambient ISM (T(J SM ). 

The different histograms in Figs. HOI and llll refer to dif- 
ferent assumptions about the available set of observations. 
The grey shaded histograms show the results obtained when 
using all the observational measurements from Fig. [9] In this 
case, all the parameters are well constrained by the observed 
ultraviolet, optical and infrared spectral energy distributions 
of the galaxies, except for the equilibrium temperature T^ c 
of warm dust in stellar birth clouds in NGC 3521. This is not 
surprising, since the total infrared emission from this galaxy 
is largely dominated by dust in the ambient ISM (as can be 
appreciated, for example, from the analogous spectrum of 



the quiescent star-forming galaxy of Fig. [S^i). We find that, 
as anticipated from Fig. [9] the transition from NGC 3521 to 
NGC 337 to Mrk 33 corresponds to a combined rise in spe- 
cific star formation rate (ips) and dust temperature (Tw° 
and Tc SM ), a larger contribution to the total infrared lu- 
minosity by warm dust (£w ) relative to cold dust (£c ot ), 
and by stellar birth clouds relative to the ambient ISM 
(ffj.), and a weakening of PAH features (^pa Vt)- These re sults 
are consistent with the finding by iDraine et alJ (|2007T ), on 
the basis of a more sophisticated physical dust model, that 
the intensity of the dust-heating starlight increases, and the 
fraction of dust mass contributed by PAHs decreases, from 
NGC 3521 to NGC 337 to Mrk 33. We emphasise the im- 
portance of finding tight constraints on all the adjustable 
model parameters in Figs. HOI and llll This confirms that the 
problem is well defined and that our model can be used to 
derive meaningful constraints on the star formation histo- 
ries and dust properties of galaxies from multi-wavelength 
observations. 

In many situations, observations may not be available 
across the full range from ultraviolet to far-infrared wave- 
lengths to constrain the physical parameters of galaxies. It 
is useful, therefore, to explore the constraints that can be 
derived on the physical parameters of Figs. HOIand llll when 
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Figure 9. Best model fits (in black) to the observed spectral energy distributions (in red) of the three galaxies NGC 3521 (top panel), 
NGC 3521 (middle panel) and Mrk 33 (bottom panel) spanning the infrared colour-colour relation of Fig. [7^.. In each panel, the blue 
line shows the unattenuated stellar population spectrum. For each observational point, the vertical error bar indicates the measurement 
error, while the horizontal error bar shows the effective width of the corresponding photometric bands. The residuals (L^° B — L™ od )/L^ bs 
are shown at the bottom of each panel. Filled squares refer to broad-band luminosties (GALEX FUV and NUV; RC3 UBV; 2MASS 
JHKs; Spitzer IRAC and MIPS; ISO, IRAS and SCUBA); open squares to the integrated Ha and H/3 luminosities. 



reducing the set of available observations. We illustrate this 
by considering three potentially common situations: 

(i) Case I (red histograms in Figs. [10] and 1110 : we re- 
lax the constraints on the Ha and H/3, /50/ISOCAM (6.75 



and 15 fim), IRAS (12, 25, 60 and 100 /jm), Spitzer /MIPS 
160 /im and SCUBA 850 fim luminosities. This leaves 
constraints from GALEX (FUV and NUV), RC3 (UBV), 
2MASS (JHKs), Spitzer/IRAC (3.6, 4.5, 5.8 and 8.0 fim) 
and Spitzer /MIPS (24 and 70 /im). 
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Figure 10. Likelihood distributions of physical quantities derived from fits to the observed ultraviolet, optical and infrared spectral 
energy distributions of NGC 3521 (top row), NGC 337 (middle row) and Mrk 33 (bottom row): fraction of the total infrared luminosity 
contributed by dust in the ambient ISM (/u); total effective V-band absorption optical depth of the dust (fy); fraction of the total 
\/-band absorption optical depth of the dust contributed by the ambient ISM (/x); effective V-band absorption optical depth of the dust 
in the ambient ISM (fj,Ty); specific star formation rate (tps)i stellar mass (Af*); and dust mass (M d ). The different histograms correspond 
to the different sets of observational constraints indicated at the top of the figure (see text of Section l3.2,2l for more detail). 




Figure 11. Likelihood distributions of physical quantities derived from fits to the observed ultraviolet, optical and infrared spectral 
energy distributions of NGC 3521 (top row), NGC 337 (middle row) and Mrk 33 (bottom row): total infrared luminosity of the dust 
(L^ ot ); global contributions (i.e. including stellar birth clouds and the ambient ISM) by PAHs (£pah)' * ne no * mid-infrared continuum 
(£m°ir) an< ^ warm dust in thermal equilibrium (C^*) to the total infrared luminosity; equilibrium temperature of warm dust in stellar 
birth clouds (T^ c ); contribution by cold dust in thermal equilibrium to the total infrared luminosity (^Q Qt ); and equilibrium temperature 
of cold dust in the ambient ISM (T<-J SM ). The different histograms correspond to the different sets of observational constraints indicated 
at the top of Fig. [TO] (see text of Section l3~X2l for more detail). 
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The most important effect of exciuding the constraints on 
Lua and Lh/3 is to iose the determination of the total effec- 
tive V-band absorption optical depth tv in Fig. 1101 This is 
because emission lines are the best tracers of dust in stellar 
birth clouds. In fact, the product fify, i.e. the part of tv aris- 
ing from dust in the ambient ISM, remains reasonably well 
constrained by ultraviolet, optical and infrared photometry 
alone. In Fig. 1111 the loss of far- infrared information redward 
of 70 jum (IRAS 100 fim, MIPS 160 jum and SCUBA 850 fim) 
has the most dramatic effect on determinations of the con- 
tribution by cold dust to the total infrared luminosity (£rj ot ) 
and on the equilibrium temperature of both cold (Tq SM ) and 
warm (T^ c ) dust. The IRAC data at wavelengths 3.6, 4.5, 
5.8 and 8.0 /im still provide valuable constraints on £paHi 
while the loss of ISO 6.75 and 15 /jm and IRAS 12 /im 
information makes determinations of f^nt and fyy more 
uncertain. The lack of far-infrared information also weakens 
the constraints on the dust mass Md in Fig. 1101 

(ii) Case II (blue histograms in Figs. [TO] and 1 1 1 p : this is 
similar to Case I, but we also relax the constraint on the 
Spitzer/MlPS 70-fim luminosity. The main effect is to sig- 
nificantly worsen estimates of the contributions to the total 
infrared luminosity by warm dust in thermal equilibrium 
(£w*) an d by the hot mid-infrared continuum (Civnk)- The 
constraints on the total infrared luminosity itself are 
also slightly worse. 

(iii) Case III (green histograms in Figs. H0land [lT]l : this is 
similar to Case II, but we also relax the constraints on the 
GALEX FUV and NUV luminosities. Not including infor- 
mation about the ultraviolet radiation from young stars has 
only a weak influence on the results of Figs. HOI and lTTI when 
optical, near-infrared and mid-infrared data out to 24 /im 
are already available. The constraints on the specific star 
formation rate (ips), the attenuation parameters (fi, tv and 
the product fity) and the relative contribution by dust in 
the ambient ISM to the total infrared luminosity (/ M ) be- 
come marginally worse. The total infrared luminosity L d ° 
and the stellar mass M* remain reasonably well constrained. 

We note that, if a galaxy has mid-infrared colours char- 
acteristic of Galactic cirrus emission, and if no other spectral 
information is available at shorter and longer wavelengths to 
constrain / M , it may be difficult to disentangle the contribu- 
tions by stellar birth clouds and the ambient ISM to the total 
infrared luminosity. Galaxies with cirrus-like mid-infrared 
emission in the SINGS sample tend to have low specific star 
formation rates and infrared spectra dominated by the emis- 
sion from ambient-ISM dust (as is the case, for example, for 
NGC 3521 in Fig. Ur). The inclusion of either far-infrared 
observations to constrain the temperature of the dust in 
thermal equilibrium or ultraviolet and optical observations 
to constrain the attenuation of starlight by dust can lift the 
ambiguity about the origin of the infrared emission (this is 
illustrated by the reasonably tight constraints obtained on 
fn for NGC 3521 in all the cases considered in Fig. [TO)) . 
We have checked that the predicted mid-infrared emission 
of the ambient ISM in these analyses is always consistent 
with the expectation that stellar ultraviolet photons (with 
A < 3500 A) are the main contributors to the excitation 
of PAH molecules and the stochastic heating of dust grains 
responsible for the hot mid-infrared continuum. This is be- 
cause stars slightly older than 10 7 yr, which have migrated 
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Figure 12. Estimates of the total infrared luminosity L^ ot de- 
rived from fits of the observed ultraviolet (GALEX FUV and 
NUV), optical (RC3 UBV) and near-infrared (2MASS JHKs) 
luminosities, plotted against estimates of the same quantity 
when including also the constraints from Spitzer (IRAC 3.6, 
4.5, 5.8 and 8.0 /im and MIPS 24, 70, and 160 /im), for 61 
galaxies with GALEX measurements in the sample studied in 
Section 13.21 The different symbols have the same meaning as 
in Fig. [7] The errors bars (in both the horizontal and verti- 
cal directions) represent the 16—84 percentile range in Lt * de- 
rived from the likelihood distributions. The bottom panel shows 
the logarithm of the difference L^ (GALEX + UBV + 2MASS) 
minus L * ot ( GALEX + UBV + 2MASS + Spitzer) divided by 
Ll ot (GALEX + UBV + 2MASS + Spitzer). 

from their birth clouds into the ambient ISM, are still bright 
in the ultraviolet. In the SINGS sample, for example, ultra- 
violet photons generated by stars older than 10 7 yr account 
for about 5 per cent of the heating of ambient-ISM dust for 
galaxies with the lowest specific star formation rates (ips)- 
This fraction reaches about 85 per cent for galaxies with 
the highest ips- At the same time, PAHs and the hot mid- 
infrared continuum are found to produce from about 1 per 
cent to about 35 per cent, respectively, of the infrared lumi- 
nosity of the ambient ISM in these galaxies. Thus, enough 
stellar ultraviolet photons are produced to account for the 
mid-infrared emission from PAHs and hot dust in the ambi- 
ent ISM of these galaxies, even if a large part are absorbed 
by cooler grains in thermal equilibrium. 

The above examples illustrate the usefulness of our 
model to interpret multi-wavelength observations of star- 
forming galaxies, and how particular observations may be 
important to constrain specific galaxy parameters. To fur- 
ther investigate the need for infrared information in the de- 
termination of -Ld"*! we compare in Fig. 1 121 estimates of this 
quantity obtained by fitting GALEX, RC3 and 2MASS data 
alone to those obtained when adding all the Spitzer mid- 
and far- infrared constraints (from IRAC and MIPS), for 61 
SINGS galaxies with GALEX measurements in our sample. 
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For most galaxies, both estimates are consistent with each 
other, although the estimates based on stellar emission alone 
are far more uncertain than those including infrared con- 
straints (as can be seen from the much larger vertical than 
horizontal error bars in Fig I12[l . There is a tendency for 
starlight-based estimates of to lie systematically under 
the more precise estimates including infrared constraints for 
the most luminous galaxies. This could arise, for example, 
if the SINGS galaxies with the largest infrared luminosities 
contained significant amounts of enshrouded stars with little 
contribution to the emergent optical and near-infrared light. 
At low L d ot , starlight-based estimates of the infrared lumi- 
nosities of blue compact dwarf galaxies in our sample are 
very uncertain in Fig 1121 We conclude that rough estimates 
of may be obtained based on ultraviolet, optical and 
near-infrared observations alone (at least for values in the 
range from a few xlO 8 to a few xlO 10 £©), but reliable esti- 
mates of this parameter require longer-wavelength infrared 
observations. We note that the total dust luminosities es- 
timated including infrared obser vations are typically within 
10 per cent of those estimated bv lDraine et alj (120071 ) on the 
basis of their more sophisticated physical dust model. 



3.2.3 Sample statistics 

We can estimate the physical parameters of all the 66 SINGS 
galaxies in our sample in the same way as exemplified in 
Figs. \W\ and [TT] above for NGC 3521, NGC 337 and Mrk 33. 
By doing so, we can explore potential correlations between 
different parameters of star-forming galaxies and infrared 
colours. In Table O we present the results of this investiga- 
tion. We list the Spearman rank correlation coefficients for 
the relations between three Spitzer infrared colours, Lf,/L 24 , 
L 2 , /L{,° an d Ll° /L u 60 , and 11 parameters derived from our 
spectral fits of SINGS galaxies: the specific star formation 
rate, ; the fraction of total infrared luminosity contributed 
by dust in the ambient ISM, / M ; the global contributions 
(i.e. including stellar birth clouds and the ambient ISM) by 
PAHs (£pah)> the hot mid-infrared continuum (£m?r) and 
warm dust in thermal equilibrium (5w*) to the total in- 
frared luminosity (equations 1191 [201 and [21) ; the equilibrium 
temperature of warm dust in stellar birth clouds, T^ c ; the 
contribution by cold dust in thermal equilibrium to the to- 
tal infrared luminosity, ^ ot (equation I22[) ; the equilibrium 
temperature of cold dust in the ambient ISM, Tq SM ; the star 
formation rate averaged over the past 100 Myr, tp = M„ips 
(see equation [34j ; the stellar mass, M»; and the ratio of dust 
mass to stellar mass, Md/M* (see equation I30p . We also in- 
dicate in Table [2] the significance levels of these correlations 
for our sample size of 66 galaxies. 

Table [5] shows that the infrared colours L% / L 2i and 
L~l°/Ll, 60 correlate well with several model quantities, such 
as Vs, /„, £pah, £w ot , T\v C , Cc ot and T C ISM , while L™/Ll° 
does not. Furthermore, it is interesting to note that the spe- 
cific star formation rate, but not the star formation rate 
itself, correlates well with Ll/L 2 V A and L™/Z/J 60 . This is il- 
lustrated in Fig. [H where we show 1%/ ' L 2 ^ and L™/Li 60 
as a function of both tps and tp for the 66 galaxies in our 
sample. We find that, as anticipated in Section 13.2.21 the 
correlation between infrared colours and specific star forma- 
tion rate arises from a drop in the relative intensity of PAHs 
and a blueshift of the peak infrared luminosity (i.e. a rise 



in the overall dust temperature) when the star formation 
activity increases. 

In Fig. 1141 we explore the relations between tps and 
other physical properties of the SINGS galaxies, which cor- 
relate well with infrared colours: / M , f;pAH, £w\ £c ot and 
Md/M* (Table [2)1 ■ The Spearman rank coefficients for these 
correlations are rs = —0.764 (indicating a 6cr significance 
level for the sample size), -0.551 (4a), 0.739 (6a), -0.647 
(5a) and 0.780 (6a), respectively. For completeness, we also 
show in Fig. [14] the relation between dust-to-gas mass ra- 
tio, Md/Mn, and ips, for the 35 g alaxies with gas ma sses 
M H = A/(Hl + H 2 ) available from lDraine et alj I (|2007h . In 
this case, the Spearman rank coefficient is rs = —0.537, in- 
dicating a 3a correlation only for the reduced sample size. 
We note that the dust masses in Fig. I14|are typically within 
50 per cent of those estimated bv lDraine et alj (|2007l ) on the 
basis of their more sophisticated physical dust model. 

The strong correlations between infrared colours and 
several physical quantities of the SINGS galaxies in Fig. [13] 
and [14] provide important insight into the link between star 
formation and ISM properties in galaxies. For example, we 
also find a la correlation between the star formation rate 
tp and the dust mass Md for the 66 SINGS galaxies in our 
full sample, and a 5a correlation between ip and the gas 
mass Mh for the 35 galaxies with H I and H2 measurements 
(not shown). To some extent, the rise in the contribution by 
stellar birth clouds (i.e. giant molecular clouds) to the in- 
frared emission from a galaxy when ^>s increases (Fig. 114a ). 
the accompanying rise in the contribution by warm dust 
(Fig. Hib). the drop in that by cold dust (Fig. ITIHI and the 
weakening of PAH features (Fig. 114b) have been anticipated 
in several previous studies (e.g..|Heloul 19861: Cesarskv et al.l 
1 19961 : ISilva et al.|[l998l : iDale et alj|200ll . 120071 ). The original- 
ity of our approach is to quantify these effects by means of a 
simple but versatile model, which allows statistical studies 
of the star formation and dust properties of large samples of 
galaxies. It is worth emphasising that our method does not 
introduce these relations a priori and that they arise from 
our consistent treatment of stellar populations and dust in 
galaxies. 

3.3 Comparison with previous models 

The model we have developed allows one to interpret the in- 
frared spectral energy distributions of galaxies consistently 
with the ultraviolet and optical emission, in terms of com- 
bined constraints on the star formation and dust proper- 
ties. At the same time, the model is versatile enough that it 
can be used to statistically derive such properties for large 
samples of observed galaxies. So far, the tools most widely 
used to extract minimum physical quantities from infrared 
spectra of large samples of obse r ved ga laxies are the spec- 
tral li braries of ICharv fc ElbaS <|200ll ) and bale fc Heloul 
(2002). As an illustration, we compare here the constraints 
derived from our model to those that would be obtained us- 
ing these libraries for the three galaxies studied in detail in 

Sect ion \3~T2\ (Fig. 191) 

ICharv fc Elba3 d 200 ll ) propose a library of template 
spectra to reproduce the observed correlations between the 
total infrared luminosity (Lir in their notation) and 

the luminosities in individual ISO and IRAS bands, for lo- 
cal galaxies. This leads them to assign a unique luminosity 
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Table 2. Correlations between three Spitzer infrared colours L^/L^ 4 , Lj 4 /-L™ and L™/Lj, eo and the median likelihood estimates of 
several physical parameters constrained using our model, for the 66 SINGS galaxies in our sample. For each combination of infrared 
colour and parameter, the first row indicates the Spearman rank correlation coefficient rj of the relation between the two quantities, 
while the second row indicates the significance level of the correlation for the sample size. 
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Figure 13. Spitzer infrared colours plotted against median-likelihood estimates of the star formation rate, for the 66 SINGS galaxies 
of the sample discussed in Section 13.21 (a) Ratio of 8-fim to 24-fim luminosity plotted against star formation rate averaged over the 
past 100 Myr, if) = M,^g (see equation [34} . (b) Ratio of 70-^tm to 160-fj.m luminosity plotted against ip. (c) Ratio of 8-/im to 24-^tm 
luminosity plotted against specific star formation rate -0s- (°0 Ratio of 70- fim to 160- /an luminosity plotted against tps- The different 
symbols have the same meaning as in Fig. [7] For each point, the vertical error bars represent the ±lcr observational uncertainties, while 
the horizontal error bars represent the 16-84 percentile range derived from the likelihood distribution of the quantity on the rr-axis. 



to a given infrared spectral shape. In Fig. 1151 we compare the 
best-fit spectral energy distributions inferred from our model 
for NGC 3521, NGC 337 a rid Mrk 33 (in black, fr om Fig.© 
with template spectra from lCharv fc El baz (2001, in green). 
Since the three galaxies have similar infrared luminosities 
(Fig. [TTJ, the y are assigned similar infrared spectra in the 
prescription of lCharv fc Elbaa (|200ll . shown as dotted green 



lines in Fig.[HJ. However, as Figs. [7]and|9]show, NGC 3521, 
NGC 337 an d Mrk 33 have differe nt infrared colours. As 
a result, the ICharv fc Elbazl (1200 ll ) template spectra that 
best fit the ISO and IRAS colours of these galaxies (plotted 
as solid green lines in Fig. I15[l correspond to widely differ- 
ent in frared luminosities. For example, the ICharv fc Elbaz 
<j2001r ) template that best fits the observed colours of the 
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Figure 14. Median-likelihood estimates of galaxy properties plotted against specific star formation rate, for the 66 SINGS galaxies of 
the sample discussed in Section 13.21 (a) Fraction of total infrared luminosity contributed by dust in the ambient ISM, / M . (b) Global 
contribution (i.e. including stellar birth clouds and the ambient ISM) by PAHs to the total infrared luminosity, £p^ H . (c) Global 
contribution by warm dust in thermal equilibrium to the total infrared luminosity £^ >t . (d) Contribution by cold dust in thermal 
equilibrium to the total infrared luminosity, §Q 0t - (e) Ratio of dust mass to stella r mass, M^/M*. (f) Ratio of dust mass to gas mass, 
M^/Mn, for the 35 galaxies with gas masses Mjj = M(Hl + H2) available from iDraine et all ll2007T ). The different symbols have the 
same meaning as in Fig. [7] The error bars represent the 16-84 percentile ranges derived from the likelihood distributions of the estimated 
parameters. 



dwarf starburst galaxy Mrk 33 corresponds to a spectral en- 
ergy distribution typical of an ultraluminous infrared galaxy 
(ULIRG, with Lir > 10 12 LfTi). This incons i stency illustrates 
the impossibility with the lCharv fc Elbad (2001) library to 
account for the intrinsic dispersion of infrared colours among 
galaxies of comparable infrared luminosity. 



bale fc Helot] (|2002h parametrise the spectral energy 
distributions of normal star-forming galaxies in terms of a 
single parameter qsed, which is inversely proportional to 
the intensity of the interstellar radiation field heating dust 
grains. Variations in cvsed can account for the observed 
range of -Fro/Fio o colours of normal star-forming g alaxies. 
iDale et all i|2007h find that the lDale fc Heloul (j2002r i models 
that best fit the MIPS colours of NGC 3521, NGC 337 and 
Mrk 33 have qsed = 3.32, 2.23 and 1.44, respectively (these 
models are shown as blue lines in Fig. I15[l . This sequence 
in osed is consistent with that in the specific star forma- 
tion rate ^>s inferred from our analysis of these galaxies (see 
Fig. I10p . As mentioned above, our model can provide greater 
insight into the star formation histories and dust properties 
of galaxies. 



3.4 Potential sources of systematic errors 

3.4-1 Star formation prior 

We have checked the sensitivity of our results to the prior 
distributions of star formation histories in the model library 
of Section 13.11 We have repeated our results using contin- 
uous star formations histories (i.e. omitting superimposed 
stochastic starbursts). We find that, in this case, the over- 
all quality of the fits to the ultraviolet, optical and infrared 
spectral energy distributions of SINGS galaxies decreases, 
and that the star formation rates of the galaxies are system- 
atically lower by up to 40 per cent. We have also t ested the 
influen ce of the initial mass function. Adopting a ISalpeterl 
(|l95Fj ) IMF instead of the Chabrier IMF would lead to stel- 
lar mass estimates about 1.4 times larger. 

3.4.2 Attenuation law 

In this study, we have assumed that the effective dust ab- 
sorption optical depth in stellar birth clouds scales with 
wavelength as rf c cx A~ 1,3 (see Section [2.1[) . Adopting in- 
stead f A BC cx A- ' 7 (i.e., the same wavelengt h dependence as 
the attenuation law in the ambient ISM; see ICharlot fc F^U 
2000) does not change significantly the overall quality of 
the fits to the broad-band spectral energy distributions of 
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Figure 15. Fits by various models to the observed infrared 
{Spitzer IRAC and MIPS, ISO, IRAS and SCUBA) spectral en- 
ergy distributions of the same three galaxies as in Fig. [9] (in red). 
In each panel, the black lin e shows the best - fit mo del from Fig. [9] 
the dotted green line is the lCharv fc Elbaj ll200ll) template spec- 
trum corresponding to the same tot al infrared luminosity as the 
black line; the solid green line is the lCharv fc Elbaj j20Qll) tem- 
plate that best fits the IRAS and ISO observ ations of the galaxy; 
and the blue line is the lDale fc H clou (200^) model that best fits 
the MIPS observations of the galaxy (see text for detail). 

the SINGS galaxies. However, for the galaxies with spec- 
troscopic observations, the Ha/H/3 ratios predicted using 
rf oc A~ 0,7 tend t o be lower than o bserved by typically 
15 per cent (see also IWild et al.ll2007l ). This has negligible 
implications for the constraints on the dust parameters /i 
and fy derived for the SINGS galaxies. 

3.4-3 Inclination 

The predictions of our model are averaged over viewing an- 
gles. In practice, however, because of the non-uniform spatial 
distribution of dust, observed fluxes at ultraviolet, optical 
and infrared wavelengths may depend on the angle under 
which a star-forming galaxy is seen. We account to some 
extent for this effect in our model when introducing the un- 
certainty <5/ M in the connection between stellar and dust 
emission (see Section [3TTJ . Also, we do not find any system- 
atic trend in the quality of spectral fits with the inclination 
of the SINGS galaxies. In the future, we plan to investigate 
the effect of orientation further by applying our model to 



larger samples of galaxies spanning wide ranges in physical 
properties and inclinations. 



3.5 Applicability of the model 

Our model is designed primarily to interpret ultraviolet, 
optical and infrared observations of large samples of (ran- 
domly oriented) galaxies in terms of effective (i.e galaxy- 
wide) physical parameters, such as star formation rate, stel- 
lar mass, dust mass, total infrared luminosity, breakdown of 
this luminosity between different dust components and be- 
tween star-forming clouds and the diffuse ISM (Section l3.2l) . 
Constraints on these quantities for statistical samples of 
galaxies at various redshifts are expected to provide use- 
ful insight into the processes that were important in the 
evolution of the galaxies we see today. 

Our model is not optimised to interpret in detail multi- 
wavelength observations of individual galaxies, for which 
the consideration of geometrical factors becomes important. 
Studies of this type generally require complex radiative 
transfer calculations for specific optical properties of dust 
grains and specific spatial distributions of stars and dust. 
We refer the reader t o the m ore s ophisticated models of , for 
example, ISilva et all l| 19981 ) and iPopescu et ail (|2000l ) for 
such purposes. Also, for detailed investigations of the ambi- 
ent physical conditions of gas and dust in starburst galaxies, 
where ori entation effect s are e xpected t o be less crucial, th e 
model of iDopita et all (|2005l . see also iGroves et all 120071 ). 
which includes radiative transfer in expanding Hll regions, 
represents a suitable alternative. 



4 SUMMARY AND CONCLUSION 

We have developed a simple but versatile model to interpret 
the mid- and far-infrared spectral energy distributions of 
galaxies consistently with the emission at ultraviolet, op- 
tic al and near-infrared wave lengths. Our model relies on 
the lBruzual fc Charlotl (|200Sl ) population synthesis code to 
compute the spectral evolution o f stellar popula t ions, and 
on the two-component model of ICharlot fc Falll (|2000T ) to 
compute the total infrared luminosity absorbed and reradi- 
ated by dust in stellar birth clouds and in the ambient ISM. 
We distribute this total infrared energy in wavelength over 
the range from 3 to 1000 /jm by considering the contribu- 
tions by four main dust components: PAH emission, mid- 
infrared continuum emission from hot dust, warm dust with 
adjustable equilibrium temperature in the range 30-60 K 
and cold dust with adjustable equilibrium temperature in 
the range 15-25 K. We keep as adjustable parameters the 
relative contributions by PAHs, the hot mid-infrared con- 
tinuum and warm dust to the infrared luminosity of stellar 
birth clouds. Cold dust resides (in an adjustabe amount) 
only in the ambient ISM, where the relative ratios of the 
other three components are fixed to the values reproduc- 
ing the observed mid-infrared cirrus emission of the Milky 
Way. We find that this minimum number of components is 
required to account for the infrared spectral energy distribu- 
tions of galaxies in wide ranges of star formation histories. 

We have generated a comprehensive library of model 
galaxy spectra by combining a library of attenuated stellar 
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population spectra (built from stochastic star formation his- 
tories and dust contents) with a library of infrared emission 
spectra. As Fig. [9]illustrates, these models provide appropri- 
ate fits to the observed ultraviolet, optical and infrared spec- 
tral energy distributions of nearby galaxies in the SINGS 
sample, for which data are available fro m GALEX, RC3, 
2MA SS, Spitzer, ISO, IRAS and SCUBA (|Kennicutt et all 
2003). We have used this model library to derive median- 
likelihood estimates of the star formation rate, stellar mass, 
dust attenuation, dust mass and relative contributions by 
different dust components to the total infrared luminosity 
of every SINGS galaxy in our sample. The accuracy of these 
estimates depends on the available spectral information. We 
find that, for example, although the total infrared luminosity 
Lj ot of a galaxy can be roughly estimated using ultraviolet, 
optical and near-infrared data alone (at least for values in 
the range from a few xlO 8 to a few xlO lo L0), reliable es- 
timates of this parameter require infrared observations at 
longer wavelengths. 

A main advantage provided by our model is the capacity 
to study the relation between different physical parameters 
of observed galaxies in a quantitative and statistically mean- 
ingful way. We find that, for example, the specific star for- 
mation rate of SINGS galaxies correlates strongly not only 
with observed infrared colours, but also with several other 
properties of the galaxies, such as the fraction of total in- 
frared luminosity contributed by dust in the ambient ISM, 
the contributions by PAHs, warm dust and cold dust to 
the total infrared luminosity and the ratio of dust mass to 
stellar mass. These correlations provide important insight 
into the link between star formation and ISM properties in 
galaxies. In particular, they allow one to quantify the re- 
lations between star formation rate, gas mass, dust mass, 
stellar mass, dust temperature and distribution of the dust 
between giant molecular clouds and the diffuse interstellar 
medium. Studies of these relations at different redshifts will 
have important implications for our understanding of the 
physical processes that dominate galaxy evolution. 

Our model should be useful for interpreting data from 
any modern galaxy survey at ultraviolet, optical and infrared 
wavelengths. It can also be used to design new observations 
by optimizing the set of observables required to constrain 
specific physical parameters of galaxies. This model is meant 
to be used by the astronomical community and we intend to 
make it publicly available. 
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